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ABS TRACT 


This thesis presents the application of eigenvalue assign- 
ment techniques to both proportional and proportional plus integral 
control of multivariable linear systems using state feedback. 

For proportional control system design, recursive design 
methods based on the eigenstructure of the open-loop system are 
considered in order to gain insight into the effects of available 
désign freedom on the performance of the closed-loop system. The 
standard eigenvalue assignment techniques developed for proportion 
control systems are modified for the design of proportional plus 
integral controllers. A new theoretical result concerning a rank 
condition for the integral control matrix is presented as a formal 
proposition and proved. 

Digital simulation studies of the eigenvalue assignment 
techniques applied to a double effect pilot plant evaporator model 
are described. The simulation results using proportional control 
schemes demonstrate that manipulation of eigenvalues only is not 
sufficient to ensure satisfactory performance of the closed-loop 
System; well-distributed closed-loop eigenvectors in the state 
Space are also essential. The simulation results also reveal that 
available design options in recursive design methods have significant 
effects on both the closed-loop dynamics and the resulting feedback 
matrices. 

Two modified eigenvalue assignment techniques for propor- 


tional plus integral control system design are applied to the 
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evaporator system. The simulation results of the evaporator system 
demonstrate that offsets in the selected state variables can be 
eliminated for any constant disturbance using proportional plus 
integral controllers designed by eigenvalue assignment techniques. 
However, the simulation results show that the transient behavior 

of the closed-loop system depends highly on the type of disturbance 


and the design method used. 
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CHAPTER ONE 


INTRODUCTION 


The dynamic behavior of many chemical processes can be repre- 
sented by a set of first-order ordinary differential equations which 
are linearized around normal operating conditions. The resulting 
state space representation of the system is in general a set of linear, 
time-invariant, multi-input, multi-output equations. It is well-known 
that the stability of such a system is exclusively determined by its 
eigenvalues and that the dynamic behavior of the system is governed 
by the modes (eigenvectors and eigenvalues) of the system [50]. 

In view of the important role of eigenvalues and eigenvectors 
in determining system dynamics, an obvious design approach is to 
synthesize a control system which ensures that the eigenvalues of the 
closed-loop system have specified values. For the class of proportional 
State feedback controllers, numerous design approaches for assigning 
closed-loop eigenvalues have been reported in the literature. Of these, 
two distinct approaches, one of which is based on the phase variable 
canonical representation of the system [1] and the other one based on 
the Jordan canonical representation (modal state representation) [41], 
42] have received the most attention. Since steady state errors (off- 
sets) due to unmeasured sustained disturbances are usually undesirable 
in process control, a natural extension of the eigenvalue assignment 
approach is the design of proportional plus integral controllers [7, 

HB rey ie 7 Oa SB 9 design of incomplete state feedback or output feed- 


back control systems is currently under active research since not all 
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the states are accessible for measurement in many real processes [3, 6, 
14]. However, the scope of this thesis is restricted to state feedback 
control systems and design approaches which are based on the modal repre- 
sentation. In addition to conceptual simplicity in design and imple- 
mentation, this approach offers the significant advantage of giving 
more insight into system dynamics and providing more design freedom. 

An important result of linear system theory [49] is that if 
a linear time-invariant system is state controllable, then state feed- 
back control can be used to arbitrarily assign all of the eigenvalues 
of the closed-loop system. Despite the many design techniques avail- 
able in the literature, reports of simulation or experimental studies 
of eigenvalue assignment techniques applied to process control problems 
are few in number. In this thesis, eigenvalue assignment techniques 
which employ state feedback are applied to a pilot plant double 


effect evaporator in digital simulation studies. 


1.1 OBJECTIVES OF THE STUDY 

The objectives of this study can be divided into two parts. 
The first part is intended to evaluate several eigenvalue assignment 
techniques for designing proportional controllers using state feed- 
back. The methods of eieanvalae assignment considered in this study 
are variations of the approach of Simon [41, 42] and will be referred 
to as "Eigenvalue Assignment via Modal Analysis" throughout this thesis. 
This method is distinguished from other methods by having a feedback 
gain matrix of the form, K = iba ge where 85 and. £, are 


column vectors of appropriate dimensions. Superscript T denotes 


the transpose of a matrix and p is a positive integer which is 
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selected a priori by the designer. According to Simon's approach, 
for an arbitrary set of vectors {g,} which preserves certain 
properties of the system, there exists a set of vectors {f,} which 
enables the desired eigenvalue assignment. It is important to note 
that the feedback controller which assigns the desired eigenvalues 
to the closed-loop system is in general not unique. 

It is the purpose of this first part to investigate the 
effects of the design freedom available in the choice of the set of 
arbitrary vectors {g,} on the performance of the closed-loop system. 
More specifically, the following factors are evaluated via the 
simulation study: 

- The effects of the closed-loop system matrix on the dynamic 
behavior of the system, e.g., if the eigenvalues of the 
closed-loop systems are the same, to what extent will the 
dynamic response of the closed-loop system change for 
different feedback control matrices? 

= The effects of the desired closed-loop eigenvalues on the 
dynamic behavior of the system. 

- The effects of the magnitude of feedback gains on the 
dynamic response of the system. 

- The effects of the set of arbitrary vectors {g,} on the 
dynamic response of the system and the resulting feedback 
gain matrix. 

= Comparison of control laws derived by eigenvalue assignment 
and optimal control theory when both have the same closed- 
loop eigenvalues. 


A further objective is to establish practical guidelines for choosing 
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the set of vectors ig,} in such a way that design objectives such 
as low feedback gains, and satisfactory transient and steady state 
behavior of the system are satisfied. 

In the second part of this study, existing eigenvalue assign- 
ment techniques are extended to the design of multivariable proportional 
plus integral controllers which eliminate steady state errors in some 
of the states. The condition that the integral controller matrix 
of the system must have full rank for arbitrary eigenvalue assignment 
is presented and proved in the form of a proposition. Based on this 
proposition, two practical algorithms are derived for designing 
proportional plus integral controllers. The performance of the 
resulting proportional plus integral feedback controllers is then 
compared with that of optimal proportional plus integral feedback 
controllers. 

In both of these studies, the control laws are evaluated 


in simulation studies by applying them to a fifth-order state space 


model of pilot plant double effect evaporator. 


1.2 STRUCTURE OF THE THESIS 

This thesis consists of five chapters including this 
introductory chapter. In Chapter Two, a literature survey of eigen- 
value assignment techniques for the design of proportional controllers 
is presented. Basic theorems and equations relevant to this study 
are also included in this chapter. 

Chapter Three deals with the extension of eigenvalue 
assignment techniques to the design of proportional plus integral 


controllers. A necessary and sufficient condition for the realization 
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of proportional plus integral controllers which assign arbitrary 
eigenvalues to the system is presented in the form of a proposition. 
Two different algorithms based on this proposition are also included 
in this chapter together with their respective modifications to dis- 
crete-time systems. 

Chapter Four describes the pilot plant double effect evapo- 
rator in the Department of Chemical Engineering and presents the 
results of simulation studies for proportional and proportional plus 
integral controllers. 

Finally, Chapter Five summarizes the major results and 


conclusions of this thesis. 
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CHAPTER TWO 


MULTIVARIBLE PROPORTIONAL CONTROLLER DESIGN 


2.1 INTRODUCTION 

Since Rosenbrock [40] first introduced modal control in 1962 
as a possible design aid in the control of large chemical plants, 
considerable attention has been paid to control system design by 
assigning eigenvalues to the closed-loop system. Eigenvalue Assign- 
ment via Modal Analysis, mainly due to Simon [41, 42], is a powerful 
method in the sense that considerable design freedom exists and 
insight can be gained into the effect of moving eigenvalues on the 
system dynamics. The design freedom enables the feedback controller 
to be designed in either a single calculation (i.e., "the simultaneous 
approach") or in a recursive manner (i.e., "the recursive approach"). 

A brief survey of Eigenvalue Assignment via Modal Analysis 
is presented in Section 2.2 including both theory and Selec as lento to 
process control systems. Since eigenvalue assignability is closely 
related to the concept of state controllability, the development of 
the relationship is also included in this section. Other eigenvalue 
assignment techniques are briefly reviewed, for purposes of comparison. 

In Section 2.3, basic equations and theory relevant to the 
development of this work are presented for both the simultaneous 
and recursive design approaches. However, major emphasis is put on 
the recursive design in view of the fact that more insight into the 


dynamic behavior of a system can be obtained by moving the eigenvalues 


one by one. 
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2.2 LITERATURE SURVEY 

The conventional proportional feedback control of a single- 
loop first-order system can be interpreted as a special case of eigen- 
value assignment since feedback control tends to decrease the time 
constant and hence shift the closed-loop eigenvalue. Actually, if 
a scalar is considered as an euetdincus ional vector, this simple 
control scheme is an example of "ideal modal control" in the sense of 
Rosenbrock. However, eigenvalue assignment techniques are also 
attractive for the design of multivariable control systems since the 
stability of the system is guaranteed by assigning desired eigenvalues 
to the system. 

Consider a time-invariant linear state space model and 


proportional feedback control system described by 
x(t) = A x(t) + B u(t) (= d/dt) (2-1) 
u(tiie= KB x(t) . (2-2) 


where x(t) isan n x-1 state vector, u(t) tsan om-x 1 icontrod 
vector and matrices A, B and KERB are real, constant matrices with 
appropriate dimensions. 

Ideal modal control was proposed by Rosenbrock [40] and may 
be defined as an eigenvalue assignment technique in which a certain 
number of eigenvalues of the open-loop system are changed while the 
eigenvectors and the other eigenvalues remain unaltered. To achieve 
ideal modal control, the control matrix B must consist of the eigen- 


vectors associated with eigenvalues to be changed [40]. However, 


in many applications, matrix B is determined by the process model 
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and cannot be easily altered. Thus the degree to which ideal modal 
control can be achieved depends on how well matrix B approximates 
the matrix of eigenvectors. Furthermore, there is no justification 
for preserving the open-loop eigenvectors unless these eigenvectors 
have favorable properties to begin with. However, systematic methods 
for synthesizing desired closed-loop eigenvectors are not available 
in the literature. 

Ellis and White [8] modified Rosenbrock's modal control 
and introduced an eigenvalue assignment technique based on a modal 
analysis of the system. They abandoned the idea of keeping open- 
loop eigenvectors unchanged and were only concerned with shifting a 
selected eigenvalue to a desired location without affecting the other 
eigenvalues. However, their efforts were mainly concerned with 
shifting a single eigenvalue in single-input systems. In an analog 
computer simulation, they applied their technique to steam pressure 
control of an oil-fired boiler described by an eight-state equation 
and achieved better control in comparison with the best conventional 
three-term controller. 

A fundamental question in eigenvalue assignment is under 
what conditions is it possible to assign an arbitrary set of eigen- 
values to the system of Equation (2-1) by the control law of 
Equation (2-2). Considerable research has been carried out to 
answer this important question and a historical survey on this 
subject may be found in the paper by Willems and Mitter [46]. The 
fundamental and important result is that the eigenvalue assignability 


of a system is identical to the state controllability of that system. 
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This property is believed to have been known for a long time for 
single-input systems but for multi-input systems, a complete proof 
may be attributed to the independent work of Wonham [49] and Simon [41]. 


This important property may be cited in the form of a theorem [49]. 


“the pair, “(AS-8)” is controllable 1f and only if’, 
for every choice of the set A, there is a matrix 
€ such that A+BC has A for its set of eigen- 


values." 


Although the final statement of the theorem is identical in Wonham's 
and Simon's work, they employed different approaches in the proof of 
the above theorem and these different approaches provide the basis 
of the two broad classes of eigenvalue assignment techniques that 


were mentioned in Chapter One. 


2.2.1 Use of the Phase Variable Representation 

Wonham's proof of the necessity part of the theorem was 
based on the linear transformation of the original system into the 
generalized phase variable canonical form. The transformed system 
matrix is in pseudo lower triangular form with blocks of phase 
variable canonical matrices (companion matrix forms) along the 
main diagonal. The generalized phase variable canonical form was 
used by Anderson and Luenberger [1] in their method of eigenvalue 
assignment. Wonham's major contribution can be interpreted as using 
the concept of cyclic subspaces to establish the existence of n 
linearly independent basis vectors which enable the system to be 
transformed into a generalized phase variable canonical form, if and 


a 


only if the system is completely state controllable. 
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The phase variable canonical form has been widely used in 
shifting open-loop eigenvalues in single-input systems due to the fact 
that the characteristic polynomial of the system can be directly found 
from this form. A major part of the computation involves transforming 
the original system into the phase variable canonical form 18, 19]. 
Anderson and Luenberger's approach [1] may be viewed as a generalization 
of the single-input case to multi-input systems. However, in contrast 
to the single-input case, the representation of a multi-input system 
in the generalized phase variable canonical form is not unique, and 
consequently different control laws result depending on the canonical 
forms used. The non-uniqueness of the phase suai canonical form 
in multi-input systems provides design freedom for objectives such 
as regulation of the feedback gain elements [34] and assigning complex 
eigenvalues to the closed-loop system [36]. An extension of this 
method to systems which are partially controllable was also reported 
[44]. The advantage of this approach is the capability of handling 
complex and/or repeated eigenvalues but the approach suffers from 
the lack of information about eigenvectors. This information might be 


valuable in process control problems. 


2.2.2 Use of the Jordan Canonical Form 

In contrast with Wonham's approach, Simon's proof of the 
above theorem was based on a Jordan canonical (modal state) repre- 
sentation of the system and provides direct insight into eigenvalue 
assignability. His proof of the necessity part of the theorem is 
also constructive. It first makes sure that the system has distinct 


eigenvalues and then these eigenvalues are moved to the desired 
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locations. Assuming that the eigenvalues of matrix A which are to 
be shifted are distinct, the postulated form of the feedback matrix 


is the sum of p dyadic products: 
Pp 
T 
slupa gect! (2-3) 


where p is the number of eigenvalues to be moved, 85 is an mxl 
vector, f. is an nx 1 vector. The structure of the gain matrix 
in Equation (2-3) provides a powerful tool in designing control 
system in such a way that eigenvalues are moved in a single step 

(the simultaneous approach), or recursively (the recursive approach) 


by specifying the sets of vectors {g, and if 5} independently. 


Simultaneous and Recursive Approaches 

In the simultaneous design method, each vector 85 is 
assumed to be a constant multiple of a predefined vector 8 in 
order to avoid having to solve a set of non-linear algebraic 
equations. This procedure has a net effect of converting a multi- 
input system into an equivalent single-input system. The vector 
g, may be chosen arbitrarily subject only to preserving the 
controllability of the original system. Although Power [35] 
suggested a quite lengthy and complicated method of choosing the 
vector 85° no systematic way of selecting a satisfactory & 
is available. Gould et al [11] extended this simultaneous design 
method to give an explicit gain formula for systems with repeated 
eigenvalues. Retallack and MacFarlane [39] presented an alternative 


approach which directly gives a compact form of the control law, 
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using the property of the return-difference determinant of a system. 
Other papers which adopted similar design approaches are also avail- 
ablen)95i243 -275429] é 

The form of the control law given in Equation (2-3) makes 
it possible to shift eigenvalues recursively in such a way that the 
p eigenvalues to be altered are divided into r (<p) groups, and 
each group of eigenvalues is moved to a desired location using the 
same number of g,'s and f,/s as the number of eigenvalues in that 
group. At each stage of the recursive design, the same procedure as 
was used in the simultaneous design is employed, and the intermediate 
closed-loop matrix and eigenvectors are calculated in order for use 
in the next stage. The design is completed by adding the individual 
control matrices, gift, that are calculated at each stage. 

If only one eigenvalue is changed at a time, the eigen- 
vectors of the intermediate closed-loop matrix can be updated from 
the previous eacdes eigenvectors by simple calculations. Ellis and 
White's method for eigenvalue assignment in single-input systems [8] 
may be regarded as a special case of this recursive design. Porter 
and Micklethwaite [27, 30] also reported a special case of this 
recursive design, derived independently from Simon, in which only 
one eigenvalue is moved to a desired location in each step and each 
eigenvalue is paired with a single control element at each step 
of the recursion. The pairing of eigenvalues and control variables 


can be achieved by specifying vector 85 to have only one non-zero 


element. 


In addition to the number of recursive steps involved, the 


basic difference between the simultaneous and recursive design methods 
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stems from the fact that linearly independent vector g, are used in 
the recursive design while only one independent vector &5 is used 
in the simultaneous design. Recursive design is a powerful approach 
due to its flexibility and the extra. insight into the effect of 
feedback on the system dynamics which is gained by changing small 
numbers of eigenvalues at a time. Furthermore, the following 
corollaries provide the necessity for and the basis of, recursive 
design in handling a derogatory system, in which more than one 
independent eigenvector is associated with a distinct eigenvalue. 

By contrast, simultaneous design cannot change all of the repeated 
eigenvalues in this situation due to the unity rank of the resulting 


controller matrix. 


"Corollary’3.6 [42] “lf the input is a scalar quantity , 
i.e., m= 1, then a necessary condition for the 

system to be completely state controllable (completely 
state observable) is that no two Jordan blocks contain 


the same mode. 


Corollary .3.14 [42]. Livthe pair (A, B) is completely 


state controllable and K isuanw (mux pn) o math. x, 


FB 


then the pair. (A. B 


lA 


. B) is completely state 


controllable." 


The approach of Simon and other workers for assigning 
eigenvalues, using the Jordan canonical form of a system, was 
termed "Eigenvalue Assignment via Modal Analysis" because the 


eigenvalues are shifted by analyzing the modal structure (eigen- 
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values and eigenvectors) of the open-loop system. Applications of 
modal analysis to other types of control problems are also found in 
the literature [5, 13, 23]. Levy et al [23] used modal analysis to 
provide insight into the dynamic response of a binary distillation 
column and to examine the validity of various column models. In 
their control system design, Davison and Chadha [5] and Howarth 
et al [13] selected the controls and states to be measured using 
modal information. Other than Ellis and White's application in the 
pressure control of a boiler, already mentioned, applications of 
eigenvalue assignment for the purpose of controlling a realistic 
eaten, using state feedback, are almost unavailable either in 
Simulation or experimental studies. This is believed to be mainly 
due to the difficulty in getting all the state measurements. How- 
ever, in the shift of a single eigenvalue using modal analysis, an 
interesting interpretation of the resulting gain-eigenvalue behavior 
of the closed-loop system obtained by a root-locus analysis was 
reported by Howarth et al [12]. 

For the class of incomplete state feedback control, where 
only very small numbers of controls and states are available, modifi- 
cations of Rosenbrock's modal control were successfully applied to 


change a small number of eigenvalues [4, 5, 13]. 


Ze) “AHEORY 


As in the case of most controller design methods, the theory 
of eigenvalue assignment techniques has been developed using a linear 
continuous-time model of a system. However, it is sometimes more 


convenient to represent a continuous system by a discretized model 
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based on a suitable sampling time, especially for purposes of digital- 
computer control. For this reason, the theory of Eigenvalue Assign- 
ment via Modal Analysis is presented for both continuous and discrete 


systems, 


2.3.1 Continuous Systems 


The dynamic behavior of many processes can be adequately 
approximated by the following linear time-invariant state space 


model, 
x(t) = A x(t) + B u(t) + D d(t) (2-4) 


where d(t) is a. q x J -di§ttirbanee vector, .D is a. constant n*x*q 
matrix and the other symbols are defined below Equation (2-1). 
The right eigenvector We and left eigenvector V5 


; , nen , ; : 
associated with the i eigenvalue ds of matrix A are defined 


by, 
Se eae nies 
Te aa T 
v, ATA, 2 : (2 6) 


Proper normalization of We and v5 gives the following property, 
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ae en ree 0) 


To investigate the effects of feedback control on the eigen- 
values of matrix A, a linear transformation, XO) as Wy(t), is 
performed on the system of Equation (2-4). Then the resulting 
expression is given in terms of a new state vector y(t), (assuming 


for convenience, dGl)e= 0), by 


WiC eye) teu tc) (2-9) 


with 


ee 
Hl 


i (2-10) 


lito 


From Equation (2-9), it is clear that the new state variables of the 
system are decoupled from each other and affected only by control 
vector ult). Matrix "H, "which ts Pete Sets called "the mode 
controllability matrix" of the system, shows the effect of the 
controls on the eigenvalues of the system. It is important to observe 
that if all the elements of the fie row of matrix H are zero, the 
aay eigenvalue of the system, hes cannot be changed by any control 
eeneme eine the oe mode (eigenvalue) is not controllable [41, 42]. 

In the absence of control, the dynamic response of the system 
of Equation (2-4) can be expressed in terms of eigenvectors and 
eigenvalues by [50], 

n rt t A, (t-T) 


n 
x(t) = ) <v, ,x(0)>e “w, + ) f <v,,Dd(t)>e” dtw. . (2~11) 
= ; -i’= ie |. <1 33 =i 
i=l i=l 6 


Equation (2-11) indicates that the stability of the system depends on 
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the eigenvalues, but the shape of transient response is also closely 

related to the orientation of the eigenvectors in the state space. 

Hence, if the control is mainly concerned with the regulation of some 

states, it is desirable to ite Cad control over both eigenvalues 

and eigenvectors in such a way that the combined effects of the eigen- 

values and eigenvectors on these particular states are minimized. 
After applying the control law given by Equation (2-2) to 


the system, the closed-loop system matrix C that results, is 
Copharin’ Keates (2-12) 


The eigenvalues of matrix C are determined by solving the character- 
istic equation det [AL | = Cilv=204, ,owhere a denotes an identity 
matrix of dimension n. If the first p eigenvalues of matrix A 
are assumed to be controllable and to be shifted to desired locations, 
P; (iG=" 172.2. .,D), suthens tue characteristic polynomial,“ £(x) >. of 


matrix C can be written as 


£Qdon detelilghe GI (2-13) 
' ey, ) 
= Il (A-p. II A-A_,. - 
i=1 a en] Pri 


However, from the assumed form of the feedback matrix given by 
Equation (2-3), the characteristic polynomial, f(A), becomes, 


: P 
£(A), =Pdet [AEP SA‘ ="Bp" )) « grpelion « (2-14) 
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Equations (2-13) and (2-14), the desired feedback matrix can be realized 


by solving the resulting non-linear equations for the g,'s [30).74L): 


2.3.2 Simultaneous Design 
In order to avoid solving non-linear algebraic equations, 
which results when one equates right hand sides of Equations (2-13) 


and (2-14) [41, 42], each vector 85 is assumed to be of the form of, 
So 4="67¢ (2424082 22950) (2-15) 


where g, is a specified constant vector of dimension m and 6. 
is a constant multiplier to be determined. Then the characteristic 


polynomial, f(A), can be written as 


f(A) =}deufAla 4 A - Bg 8 Vv] (2-16) 


: : ath 
wnere 6 isa “p <,l- “vector, y is an nx p matrix and the i 


element of 6 and the wig column of YS are denoted by 65 and 
Vi» respectively. Applying the Schur fomula [10] and using the 
relation given by Equation (2-8), the right hand side of Equation (2-16) 


can be factored into two determinants, hence 


£(A) = det (AZ -A] det[1-8"(L A) a] (2207) 
where 
A = diagla, sry 2-545] 
@= EEN OaN = VBE, 
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T -1 ree 
det[1-6 QI -A) “oJ =1- ) —= > . (2-18) 
ie oe ae i=l i 


Suppose that oF are specified such that the following relation holds 


Be OO. Pp (0-0) 


alae hi 
Wit beh Gacgey, ical O-A,) 


(2-19) 
ate) i=1 


Then from Equations (2-17) and (2-18), the characteristic polynomial 
£(A) , becomes identical to Equation (2-13), thus assigning desired 
eigenvalues to the closed-loop system. If a. #0, the constant 
coefficient Fi can be determined from Equation (2-19) and is given 


by 


SS. be ania Geers, wae (2-20) 


j=l 
see! 
It is observed from Equation (2-20) that the underlying assumptions 
made in the derivation are that re # °, Ci = 2 oe) cand 


that g is chosen such that vector a has non-zero elements. 


2.3.3 Recursive Design 

As mentioned in the previous section, the form of the feed- 
back matrix given by Equation (2-3) suggests the possibility of shift- 
ing eigenvalues recursively. Although, an arbitrary number of eigen- 
values can be shifted at each step of the recursion, only the case 
where one eigenvalue is moved to a desired location at each step, is 


dealt with here. From Equations (2-12) and (2-3), the closed-loop 
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can be written as 


(2-21) 


eigenvalues of matrix A, de (GR BRAN ego 9 Re 
are assumed to be shifted to new locations, P5 Pa LADS Ree pare oD Ee 
The intermediate closed-loop matrix, Ci after the yth 


recursion, is defined by 


and 


Clearly, Equation (2-22) 
the ee recursive step 
step. Assuming that the 
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prespecifying vector 85 
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shows that the closed-loop matrix c. after 
becomes the open-loop matrix for the ar 


ath d F : 
i eigenvalue dy of matrix A is moved 


step, the desired objective can be achieved by 


arbitrarily and setting 


fe de SD (2-24) 
=i i -i 
, C121) ; 
where 65 is a constant to be determined and VY; is defined by 
the following relation, 
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= -2 
¥i-1 Vi Ae (2-25) 
; Gila 3 : 
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.th : : 
the i step. By applying Equation (2-20), 65 can be calculated: 


(9 ;-A,) 
65 = a ay i aah Cd: Stag? 20. .p), (2-26) 
gee. Be. 


Observe that each vector 85 should be chosen such that 
<*> in Equation (2-26) is not equal to zero and this is always possible 
since the first p-. eigenvalues of matrix A--are controllable. The 
rank of the resulting feedback matrix is less than or equal to no, 
depending on the number of independent vectors, 85° By contrast, 
in the simultaneous design, the rank of the feedback control matrix 
is always unity. It is interesting to observe that the elements of 
vectors, 85> can be interpreted as the ratio of control efforts 
employed in shifting the ees eigenvalue [41]. If all the elements 
of vector g, are non-zero, all controls are used simultaneously and 
if only one element is specified to be non-zero, then a single control 
is used individually to shift each eigenvalue | A;° Actually, Porter's 
sequential eigenvalue assignment algorithm [27, 30] is an example of 


the latter case. 


20324 Denaeerace of the Updated Right and Left Eigenvectors 

In the recursive design method, it is ees tae to calculate 
left eigenvectors at each recursive step, but if a small number of 
eigenvalues are shifted at each step, right and left eigenvectors of 
matrix Cc. can be easily updated from those of C. 4 [41,42]. Again, 
only the case, where one eigenvalue is changed at a time, is described 
here. Although it is not absolutely necessary, the eigenvalues of 


matrix C, 1 are assumed to be distinct for convenience, and vectors 
=l1- 
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W and denote the right and left eigenvectors of 


(i-1) (i-1) 
i Vi 


Sit 
: : _th 
associated with the i eigenvalue. Let the set of eigenvalues of 


matrix Cc. after the AG recursive step, be Oso Ged RP aeeca ans 30 la 


‘ ‘1 
then ,4) = ee ) (if j # i). From the definitions of an eigen- 


vector and matrix Ci it is easily verified that, 
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By the assumption of distinct eigenvalues of C, the set of vectors 
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{ aad (j = 1,2,...,n) spans the n-dimensional vector space, and 


hence it is possible to write, 


esl oo (2-28) 


n 
wa = oy qj le 


From the definition of an eigenvector and using the property given 


by Equation (2-7), it can be shown that [41] 
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using the property given by Equation (2-7), it can be readily verified 
that 


Cyt 


j Ae CPE) sd) (2-31) 
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2.3.5 Controller Gain Minimization 

In many practical problems, it is desirable to put constraints 
on the magnitude of the controller gain elements because of limitations 
due to physical control variables, measurement noise, etc. In both 
simultaneous and recursive design methods, design freedom is provided 
by an arbitrary vector g, or a set of arbitrary vectors ig, } with- 
out affecting the eigenvalues of the final closed-loop matrix. If 
regulating the magnitude of the gain elements is a prime concern, these 
vectors can be chosen so as to minimize the magnitude of the largest 
gain element. In general, this will require the solution of a non- 
linear optimization problem. However, if only one eigenvalue is 
shifted, Simon [41, 42] showed that the largest gain element can be 
minimized by selecting vector o to be 
go = sign LS (hasta 25 lew) (2.32) 

J 

where g, is the es element of vector Bo is the left eigen- 
VECCOr, il peke with the eigenvalue to be shifted and 2 is the 
hl column of matrix B. Simon derived the relation of Equation (2-32) 
by defining "a measure of controllability" for the eigenvalue to be 


changed and proving that 8B chosen by Equation (2-32) maximizes 


this "measure of controllability". 
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Although Equation (2-32) is only valid when a single eigen- 
value is shifted, this policy can be used in a recursive design to 
preserve the controllability of the system and to prevent excessively 


high gain magnitudes. 


2.3.6 Discrete Systems 

In many cases, the actual implementation of a control scheme 
can be more conveniently achieved on a discrete-time basis rather 
than on a continuous-time basis, especially in modern digital computer 
control systems. In such a case, the continuous model given by 


Equation (2-4) can be discretized into an equivalent discrete model 


[25], 
x[(k+1)T] = 9(T)x(kT) + A(T)u(kT) + 9(T)d(kT) (2-33) 
where 
AT 
OCT) s ce 
TAPAT 
A(T) = (f e  dr)B 
Oo 
T At 
e(T)aac) ey Pd2)D 


and T is the sampling time of the system. 

It can be easily shown that if the original continuous 
system is completely state controllable, then the discrete system 
given by Equation (2-33) is also completely state controllable [46]. 
Furthermore, the right and left eigenvectors are the same as the 


2 ath 2 
corresponding quantities of the continuous system and the i eigen- 
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value of the discrete system is e - If the control is based on 


discrete control law of the form of 
u(kT) = Kop x(kT) (2-34) 


then the design techniques described in previous subsections can be 
applied without any modification. However, to insure the stability 
of the resulting closed-loop system, the eigenvalues of this discrete 


system should be shifted inside of the unit circle in the complex plane. 
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CHAPTER THREE 


MULTIVARIABLE PROPORTIONAL PLUS INTEGRAL CONTROLLER DESIGN 


3.1 INTRODUCTION 

The presence of sustained external disturbances is a common 
occurence in many practical control problems. It is well known that 
proportional control alone cannot eliminate undesirable steady-state 
errors (offsets) arising from sustained disturbances. Consequently, 
it is common practice to use integral feedback control in conjunction 
with proportional control for this purpose. 

It is the objective of this chapter to consider the design 
of a class of multivariable proportional Mine ineeer el controllers 
(which will be referred to as PI controllers by applying eigenvalue 
assignment techniques. A survey of pene texant literature is presented 
in Section 3.2 with emphasis on the application of eigenvalue assign- 
ment techniques. The relationship between eigenvalue assignability 
and the state controllability is also reviewed in this section. In 
Section 3.3, a necessary and sufficient condition for realizing a 
PI control law which enables arbitrary assignment of all system eigen- 
values is presented and proved. Two practical algorithms for assigning 
arbitrary closed-loop eigenvalues are also summarized in this section. 
In Section 3.3.2, these algorithms for continuous systems are modified 


for application to discrete systems. 


3.2 LITERATURE SURVEY 


For the class of state feedback control system, Johnson [15- 


17] was one of the first investigators who considered adding integral 
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control to proportional control to regulate states. He first considered 
the case of constant disturbances [15, 16] and excluded disturbance 
term from the state equation by adding an auxiliary state variable 
which represents the combined effect of control and disturbance terms. 
An essential assumption in this approach is that the range of the 
coefficient matrix, D, of the disturbance vector is contained in 
the range of the coefficient matrix, Bs, jot ethexcentrol:veetons (Lé 
the state equation satisfies this condition, he showed that offsets 
can be eliminated in n state variables by applying optimal control 
theory to the augmented system equation. In a later paper [17], 
Johnson extended his previous approach to the class of arbitrary 
vector disturbances, which satisfy a linear differential equation, 
and described controller design methods using both optimal control 
and stabilization (eigenvalue assignment) theory. Johnson's design 
approach using eigenvalue assignment technique guarantees the regu- 
lation of all state variables in the presence of a broader class of 
unmeasurable disturbances. However, in order for the essential 
assumption employed in his developments to be applicable to an arbi- 
trary system, it is required that the number of controls be equal to 
the number of states. In many applications, this condition is dif- 
ficult to meet because available controls are usually limited in 
number. 

A different approach to this regulation problem using 
optimal control theory is given by Porter [25] and Newell and Fisher 
[21]. Their approach is different from Johnson's in that only as many 
integrated states, as preserve the controllability of the augmented 


system, are fed back, hence eliminating the range condition imposed 
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on the coefficient matrices of the control and disturbance vectors. 
However, this approach will eliminate offsets in only r_ states 
where r<m and is valid for impulse or constant disturbances, 
Newell experimentally applied his approach to the double effect evap- 
orator in the Department of Chemical Engineering and observed improve- 
ment over a proportional control scheme. 

Before, proceeding further, it is convenient to define the 
problem mathematically. Consider the system described by Equation (2-4); 
ifaait-is, desired. to,eliminate.offsets) in, »).(<.n).. state, variables, 
then a simple treatment is possible by defining a rx 1 _ vector, 


z(t), by 


x(t) (3-1) 


where Z. is a r xn matrix which consists of r appropriate rows 
of the nxn _ identity matrix Le Combining Equations (2-4) and 


(3-1) gives the augmented system: 


> 


k(t) = A x(t) +B u(t) + D d(t) (3-2) 
where 
, Oe 
eCoye= 3 “ZCt) C3=3) 
: Ae | 
A = | -=--1---- (3-4) 
x totaled gineea! 
‘acd yoo 
(3-5) 
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(3-6) 


lit > 
1) 
t 
to) ie} 


and the O's are null matrices with appropriate dimensions. 

The effect of augmenting the state vector is the introduction 
of r repeated zero eigenvalues to the system in Equation (3-2). 
These r zero eigenvalues should be moved to the left half of the 
complex plane in order to stabilize the system. The problem of regu- 
lation of r states is now converted to assigning at least r eigen- 
values of the (ntr) x (ntr) system matrix, A. In view of the 
relation between state controllability and eigenvalue assignability 
of a system [41, 49], it is necessary to examine the controllability 
of the augmented system. The conditions for the controllability of 
the pair (A, B) have been extensively investigated and necessary and 


sufficient conditions in terms of matrices 


A, B and Z- were reported 
by Porter and Power [31, 32, 38] and Davison and Smith [7, 43]. In 
their first paper [31], Porter and Power showed that r cannot elteca, 
m, the maximum rank of acetne B and that a necessary condition for 
(A, B) to be a controllable pair is that the pair (A, B) must be 
controllable. They later prove that necessary and sufficient conditions 
for the controllability of the augmented system are that (i) the pair 


=] 


(A, B) is controllable and (ii) the rank of matrix ZC +BK) Bn sks 


r, where matrix K must be chosen to ensure the invertability of 
matrix A [32, 38]. 

Davison and Smith have investigated the controllability of 
the augmented system with extension to the general case where matrix 


T is arbitrary. They presented the same conditions given by Porter 
wis 
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and Power but in a more general form, and proved that if all the zero 

and positive eigenvalues of matrix er are moved to the left half of 

the complex plane, offsets in r states are eliminated for constant 

or impulse type disturbances. However, to cope with other types of 
disturbance such as a ramp or a parabolic disturbance CEC. arepro- 
portional plus multiple-integral control system should be considered 

[7, 17]. The necessary and sufficient conditions for the controllability 
of the augmented system of Equation (3-2), which ti abike arbitrary 
eigenvalue assignment, are given by [7]: 


i) (A, B) is a controllable pair and 


Once the controllability of the augmented system of Equation 
(3-2) is ensured,then it is possible to shift an arbitrary number of 
eigenvalues of matrix A (< ntr) by designing a proportional controller, 
ao= K x, for the system in Equation (3-2). However, if the simulta- 
neous design approach described in Section 2.3.2 is applied to this 
augmented system, at most only one of the r repeated zero eigen- 
values can be shifted due to the unity rank of the resulting feedback 
matrix [33]. This restriction is an obvious result of Corollary 3.8 
cited in Section 2.2 and implies that at most only one state 
variable can be guaranteed to have zero steady state error. On the 
other hand, the recursive design approach described ineSectionn2 2313 
can be applied to shift an arbitrary number of eigenvalues of matrix 
A to the desired locations [29]. However, it should be kept in 


mind that a minimum of r_ steps of recursion are necessary to eliminate 
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offsets in r state variables and that the (ntr) dimensional eigen- 
vectors of matrix A must be updated or recalculated at each recursive 
step even when shifting the eigenvalues of matrix A. Another possible 
approach is to apply Anderson and Luenberger's method of eigenvalue 
assignment to the system of Equation (3-2). 

The special structures of matrix A, B and D make it 
possible to decouple the task of assigning eigenvalues of the 
(ntr) x (ntr) system into two successive eigenvalue assignments for 
(n x n) and (rx xr) systems. This approach was adopted by Power 
and Porter [37] in that, after first shifting the eigenvalues of the 
original system matrix A, the r zero eigenvalues of A are moved 
to desired locations simultaneously by inverting a suitably chosen 
(cr x r) non-singular matrix. A similar design technique is derived 
in Section 3.3.1 from a condition concerning the rank of the 
integral feedback matrix. 

Only a few applications of eigenvalue assignment approach 
to the design of PI controllers have been reported in the literature 
[4, 5]. Davison [4] applied an eigenvalue assignment technique to 
design the integral controller in the PI control of a boiler system 
described by 9 states and 2 controls, and a distillation column with 
11 states and 3 controls. His simulation results showed that offsets 
were eliminated in those states subjected to integral action. Another 
application provided by Davison and Chadha [5] was concerned with the 


regulation of a single state variable of a large composite chemical 


plant of 41 states and 8 controls. 
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3.3 NEW RESULTS 

The main new result is presented in Section 3.3.1 in the 
form of a proposition; two eigenvalue assignment algorithms are 
described based on this proposition. In Section 3.3.1, continuous- 
time systems are considered and the two algorithms are modified for 


application to discrete models in Section 3.3.2. 


3.3.1 Continuous Systems 

Consider the system described by Equation (3-2). The 
following proposition provides a statement of the restriction on the 
integral feedback matrix which guarantees the arbitrary assignment 
of all ntr eigenvalues. 

Proposition: Suppose that the system represented by 
Equation (3-2) is completely state controllable. Then a linear feed- 
back control law of the form of Equation (3-7) which shifts all ntr 
eigenvalues of matrix A to arbitrary non-zero values, which are 
either real numbers or complex conjugate pairs, can peuesniie es at 


and only if the rank of the matrix, K is equal to r. 


Thad 


u(t) = [Kyp | K,]x(t) G7) 
Proof: Combining Equations (3-7) and (3-2) give the closed- 
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(Necessity) It will be shown that if the rank of matrix Kt 


a 


is less than r, the closed-loop system matrix C retains at least 
one zero eigenvalue even though the system of Equation (3-2) is com- 


pletely state controllable. Matrix C can be expressed as the product 


of an (nt+tr) x (n+m) matrix and an (ntm) x (ntr) matrix as follows, 


Since the system of Equation (3-2) is completely state controllable, 
the first partitioned matrix in Equation (3-9) will have a rank of 
(ntr). However, from the second partitioned matrix in Equation (3-9), 


it is clear that if matrix Kt has a rank less than r, then the 


“a 


rank of matrix C becomes less than (n+r) because the rank of the 
second matrix is less than (ntr). This implies that at least one 
zero eigenvalue is retained. Thus the necessity part of the pro- 
position is proved. 

(Sufficiency) This part of the proof is constructive since 
a linear feedback control law of the form of Equation (3-7) is derived. 
Assume that K_ has rank r. Denote the spectrum of eigenvalues of 


8 


r_,0,...,0} and the spectrum of desired non-zero 


oan 


A by {\,; 


eigenvalues by {Py Poorrs Py sPrggsrssoPryet: 


It is convenient to write the control vector u(t) as 
aCE) = Cee) (3-10) 


where u, (t) will be used to shift the n eigenvalues of matrix A 


and u, (t) to shift the r zero eigenvalues of A. Next, u, (t) is 
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specified to be of the form 
(t) = [K} Q]x(t) (3-11) 


where K is an mxn matrix. Combining Equations (3-10), (3-11) and 


(3-2) and assuming d(t) = 0 for convenience, gives the following 


expression 
x(t) = A)x(t) + Bu, (t) (3-12) 
where 
: A+B! 9 | 
ae.,| RR2Eee SELLSG_& = 
A, = iO : (3-13) 
=r ' = 


The eigenvalues of the above system are the roots of the 


equation 
det [AL - (AtBK)] ° det [Ar] =7Q *. (3-14) 


Equation (3-14) indicates that only the n eigenvalues of matrix A 
are affected by u(t). Since the pair (A,B) is controllable, a 
proportional feedback matrix K can be designed to assign the desired 


eigenvalue spectrum, {0 sPos+++ Pits to matrix, (A+BK). 


“a OM li 
After «K) ds specified, the modal matrix V of AL has 
the form. [32]; 
‘ 

se Go aya 
Yo = | -----=~----;-++---- (3-15) 
a <1) (AtBK) eo) aL 

=r = == ¢ =r 


T : d 
where V is the modal matrix of (A+tBK) . Performing the linear 
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T 


transformation, y(t) aa KEE) Equation (3-12) is transformed into 


Jordan canonical form, 


“ Hy 
YOM tesa os ECL) ee eeeee s(t) (3-16) 
t| a 


where J is the Jordan canonical form of A+BK and 

H - Y'B (3-17) 
H = -T.(AtBK) BO (3-18) 
The control vector u, (t) is specified to be 

(t) = [2 1 K ly¥(t) (3-19) 


Or,ein terms “of x(t)), 


lz). - (3-20) 


Combining Equations (3-16) and (3-19) gives the closed-loop matrix in 
the modal domain, Cc, , 


oe Ale. (3221) 


HO > 
tl 


The (ntr) eigenvalues of c, are the roots of the characteristic 


equation, 
det [AI -J] : det [AI -HKy] SOU ys (3-22) 


From Equation (3-22), it is evident that Kr does not affect the 
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eigenvalues of J, which reduces the problem to assigning the eigen- 
values: of the’ r-* r' matrix, HK: 


Matrix K is specified to be of the form 


it 


ig 
Kpeatneg tf (3-23) 


where 85 is an arbitrary m x 1 vector subject only to the control- 
lability requirement and ft. isan or <x L syector.) Since it ‘is 
assumed that matrix is hasitrank: relia both {g,} and {f, should 


be specified so that linearly independent vector sets result. A rxr 


matrix Cc. is defined recursively by 


a ii ha 
ae. Cet eee ie 9 eorels ork tot) (3-24) 


with 


(3-25) 


ta 
i] 
"lo 


Then, the recursive design method described in Section 2.3.3 can be 


directly applied to shift the r zero eigenvalues of Ss by 


specifying 
a, ee (3-26) 
=i i-i 
(i= en ; 
where v, is defined as the j left eigenvector associated with 
ma 


(o) 


the j h eigenvalue of C.> however, from the Equation (3-25), . v, 
aan be chosen arbitrarily but the set fv nie should be linearly 


independent. 
By combining Equations (3-10), (3-11) and (3-20), the 


resulting control law can be expressed as 
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Cts o[ kek PCL) rire (3-27) 


This proposition serves as the basis for two algorithms for 
assigning eigenvalues. In both algorithms, the first step is to shift 
selected eigenvalues of A to desired locations using proportional 
feedback control. The two proposed algorithms then differ in the 
design of the integral feedback matrix, Kie which is used to shift 
the r zero eigenvalues. 

Recursive Design: This algorithm closely follows the proof 
of the sufficiency part of the proposition in that the r zero eigen- 
values are shifted to desired locations recursively. The algorithm 
consists of the following three steps. 

Step 1) Design the m <n _ proportional control matrix K 
so that matrix, (A+tBK), has the desired eigenvalues. This can be 
achieved using the algorithms mentioned in Chapter Two. 

Step 2) The integral feedback matrix, Ki -is specified 


to be of the form 
r P T 
Kane y ie! : (3-28) 


5H g,> 7 O- “and theset, {gts is linearly 


(i-1) 
Vi 2=i 


Choose 85 so that < 
independent. Since the augmented system of Equation (3-2) is completely 
state controllable, it is always possible to choose vectors 85 which 


satisfy these conditions. As shown in Section 2.3.3, 65 can be 


calculated from the equation 


ree comes AMR ea 5 Cae ee (3-29) 
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A systematic way of specifying 85 is to follow the procedure des- 
cribed in Section 2.3.5 which ensures the controllability of the 
system. 

Step 3) Complete the design by constructing the feedback 
control law given by Equation (3-27). 

This algorithm retains the computational advantage of the 


unity rank (dyadic) approach but yields a i matrix of rank. it, 


which shifts the r zero eigenvalues. The left eigenvectors, Acie 
which are required at each stage, can be updated by the method des- 
cribed in Section 2.3.4. 

Simultaneous Design: Steps 1) and 3) of this algorithm are 
identical to those of Recursive Design. Consequently, only Step 2) 
is presented. 

If the control law of Equation (3-27) is used, it follows 
from Equation (3-22) that the (n+tr) eigenvalues of the closed-loop 
system consist of the n_ roots of det [Al -(AtBK) ] = 0 andthe r 
roots of det [AI -H,K,] = 0. This last equation can be used to 
design Kie 


Step 2) Let 


Kole Dp (3-30) 


where D is a rx r matrix whose eigenvalues are the r desired 


eigenvalues, 0 bd (i = 1,2,...,r). Furthermore, assume that KY 


has the form 
T 
ee acting (3-31) 


where G is a rr non-singular matrix. From Equations (3-30) and 
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(3-31), it follows that 


Kop) 
NW 
io 


(3-32) 


Since the rxm matrix Cis.m) 5 Ha» Naserank sh saethengthe. (xr x <r) 
: 2 
matrix HAH, also has rank r. Hence, Equation (3-32) can be rearranged 


to give 
Gr = (HH) Denes (3-33) 


Combining Equations (3-31) and (3-33), matrix K. is expressed by 
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The apparent difference between two algorithms is that the 
Recursive Design requires the calculation of eigenvectors, while the 
Simultaneous Design employs matrix inversion. This distinction is 
characteristic of all eigenvalue assignment techniques in the 
literature in that either eigenvector calculation or matrix inversion 
is required. However, the choice between the above two algorithms 
should be based on other factors such as magnitude of controller gains or 
performances of resulting systems rather than upon the computational . 


difficulties involved. 


Comparison With Existing Methods 


As mentioned in the last part of Section 3.2, if the 
augmented system of Equation (3-2) satisfies the controllability con- 
ditions, eigenvalue assignment techniques developed for proportional 
control systems, such as Anderson and Luenberger's approach [1] and 


recursive design based on modal analysis [26, 41, 42], can be success- 


Cet ay fy q - 


(@£-£) ee a = 8 leon 


Yd boeeeardxs et 


«2 wtvisa) , (@é-t) bins WEE) enolyaunag 


(sE-E) 


ons elisiw «S107 osvaagte to no bie Ine leo” ak ty oe 
ab noksontasth eidt. vobkesavad xa 
sil ok pouptauiaes insmog tees suley " 


aoledevat x i3em 10, nolan Lipo kao" aaieniis 


10 entey shad ccxceaens't do sbucitoaga a foue ex03502 alte 10 


i" 


ees arta oqu: neds reioer amaseye es toe 


Geet 


=) 


Sai tb S-E: aotsoa? 


ay 


9 AW 338g, tae! a at 
. “mao Wabldestoxss00 oid. ‘eoitetaes (Se) © ont i tc 


40 


fully applied to assign all eigenvalues of the augmented system. This 
explains why investigators in the field of eigenvalue assignment have 
been mainly concerned with establishing controllability conditions for 
the augmented systems. However, the disadvantage in using those 
existing techniques is that one has to deal with (ntr) dimensional 
Systems rather than n dimensional systems. Although r is usually 
much less than n, the computational burden increases drastically as 
the number of states increases. Actually, if as many controls as 
States are available, then probably Johnson's approach [17] is the 
most useful approach since it can cope with a general class of 
disturbances and the dimension of the augmented system is always 

(n+1) if constant disturbances are assumed. 

The two design algorithms derived in this section are based 
on the strategy of decoupling the eigenvalue assignment problem for 
the (ntr) x (ntr) system into two smaller problems, namely, assign- 
ing the eigenvalues of an n x n system and of a r xr. system. This 
strategy is not new and has been used in the design of PI controllers 
by Johnson [17] and Power and Porter [37]. Its chief advantage is 
that only eigenvectors of the smaller subsystems need be calculated 
rather than the eigenvectors of the (ntr) x (ntr) system. 

The second method given in this section is similar to that 
of Power and Porter [37] except a pseudo-inverse of matrix H, is 
used in the former, while the inverse of a matrix consisted of r 
linearly independent columns of matrix Hy is used in the latter, 
for the design of the K, matrix. If m is equal to r,° these 
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two approaches will give identical results. 
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In the first method described in this section, although the 
basic design scheme is the same as used in Power and Porter's approach 
[37], xr stages of recursion are used in designing the integral con- 


troller based on Corollary 3.8 given by Simon [41, 42]. 


3.3.2 Discrete Systems 

If a discrete PI control law is preferred over a continuous 
one, the discrete controller may be derived by applying the eigen- 
value assignment techniques to a discrete equivalent of Equation (3-2). 
Different discrete representations can be derived by several approaches. 
For example, one can directly discretize Equation (3-2), following 
Standard procedures in the literature [25]; this gives the following 


equations: 


x[(kt1)T] = $(T)X(kT)+A (T) u(kT)+8(T) d (kT) (3-35) 
where 
; x (kT) 
x (kT) = ~Z(kT) ° (3-36) 


Matrices $(T), 8(T) and A(T) are defined in the same way as 
matrices $(T), Q(T) and A(T) are defined in Section 2.3.6. How- 


ever, if A is non-singular, matrices o(T), 6(T) and 4(T) can 
+ 


be expressed compactly by 


$(T) = | nae ceg-nan----- tae a, (3+37) 


The details of the derivations of Equations (3.37)-(3-39) are 
presented in the Appendix. 
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(3-38) 


(3-39) 


As an alternative strategy, the discrete representation of 


the augmented system can be derived by combining Equation (2-33) with 


the following equation which defines a rx 1 vector, 


or equivalently 


2 ((kt1)T) = TOT x[ (k+l) T] +2 (kT) 


z(kT), by 


(3-40) 


(3-41) 


Combining Equations (2-33) and (3.41) gives the augmented discrete 


version of the system, which has the same form of Equation (3-35). 


However, matrices 4(T), (T) and A(T) are now modified to 
: g(T) , Q 
sO) = oar gay f ie 
P 8(T) 
8D) = |-aE-9(ay tes 
A(T) 
= (3-44) 
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Newell and Fisher [21] used Equations (3-42) to (3-44) in their optimal 
PI controller design and experimentally proved the utility of these 
equations by applying the PI controllers to a pilot plant double effect 
evaporator. 

Another discrete representation of the augmented system can 
be found in Porter and Crossley [29]. Their alternative definition of 


z(kT) is given by 


k-1 
ZCKT) =eToe) Text) (3-45) 
izo * 
or equivalently 
z[(kt1)T] = TE x(kT)+z(kT) . | (3-46) 


Obviously, matrices 4(T), (T) and A(T) are modified and now given 


by 
, $(7)10 
g(T) = -uaiaN (3-47) 
) =r j=r 
To §(T) 
CLG ih da pata = (3-48) 
. A(T) 
NOW a Ir icity aia he (3-49) 


It is observed that in these three alternative discrete 
representations of the augmented system, the augmented system matrix 


$(T) has the same set of eigenvalues. Newell and Fisher's and 


fm ge abi ak (osteo to) * 
| séad3 to) eatiian: ot. bevers. ‘ | 
ta9326 ebduob sisig doltg 03 erefl 
f25 wereye besuemaen ex 22 noddeimesosson eter 
40 noidinbish ovttsayssis sPecT is) ‘weleeot baw. fentot at ba 


{éa-E) : (Tada ( 9 = eng: | i 7 ; 


AA L 

ows, 

t A Mv ae, 
oe ae : 

ul 


(OBE) . sh Dar “ee 


novig wor bas beiitbom sis” (Dd | bats “Oe. ont tetas 


44 


Porter and Crossley's expressions for z (kT) can be interpreted as 
difference approximations of Equation (3-1). The only difference in 
these two representations is that a backward difference equation is 
used by Newell and Fisher while Porter and Crossley used a forward 
difference equation in discretizing Equation (3-1). Obviously, the 
accuracy of these representations will depend upon the magnitude of 
the sampling time, T, used. For a sufficiently small sampling time, 
Porter and Corssley's representation will give the simplest version 

of the continuous augmented system. However, in this study, Newell's 
representation has been adopted since this approach was experimentally 


tested. 


Discrete PI Controller Design 


It is clear from Equation (3-42) that the eigenvalues of 
the augmented discrete system is the union of the eigenvalues of the 
nxn matrix $(T) and the rx r_ identity matrix I. For 
satisfactory PI control, the r nitty eigenvalues should be moved 
to inside of the unit circle in the complex plane by employing a 


discrete control law of the form: 


u(kT) = [K,,{ K,Jxr) . (3-50) 


The complete state controllability of the system given by 
Equation (3-35) is a necessary and sufficient condition for arbitrary 
eigenvalue assignment in view of the relation between state control- 
lability and eigenvalue assignability of a system [41, 42, 49]. These 
controllability conditions for the augmented discrete system of 


Equation (3-35) can be stated as: 
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i) thevpaic” [o(1). A(T)] is a controllable pair, 
1 
ii) rank | ----=---- trcacco- = ntr 
i) 


This can be proved by using the argument given by Davison and Smith [7] 
for continuous systems. 

If the augmented system of Equation (3-35) satisfies these 
two conditions, the proposition and the two algorithms described in 
Section 3.3.1 can be extended to shift as many of the ntr_ eigen- 
values as desired, to arbitrary specified positions. Due to the 
difference in the structure of the augmented system matrices for 
continuous and discrete systems, some modifications are necessary for 
two design algorithms to be applicable to the augmented discrete 
system. Only those modifications, other than simple substitution of 
o(T) and A(T) for A, B, are described below. 


For discrete systems, Equation (3-15) should be modified to 


Vee jh ee oe ee ; (3-51) 
“TT (Lt1gDty@e-L, I} j 
By performing a linear transformation on Equation (3-35), y(kT) = 
Vi x(kT), a Jordan canonical form similar to Equation (3-16) is 
obtained with the identity matrix i. in the lower diagonal block 


and Equation (3-18) becomes 


-1 
H, = -TT[e(T)+4(DK-L J ac). (3-52) 


With this modification of matrix Ha» Equation (3-29) becomes 
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Pret aes 
(i512) f 
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(3=53) 


The location of the r eigenvalues which are moved in Step 2) of the 
Simultaneous Design is determined by the roots of det [AL -1,-E,K,1 = 0. 


Hence, the necessary modification of Equation (3-32) is given by 
Re Ds (3-54) 


Then, clearly the integral controller matrix K is represented by 


E 


Equation (3-34) becomes 


a T.-1 
HB apt Del dy. (3-55) 


WA 
i} 


Finally, the closed-form of the feedback PI controller matrix is 


given explicitly by 


ms oe ; 
SUKOTE Ad Ela D eA) Koes Mh. (3-56) 


WA > 


Equation (3-50) is the discrete equivalent of Equation (3-27) for the 
continuous augmented system. 

It is emphasized that the general procedure which is des- 
cribed in Section 3.3.1 is applicable to both continuous and discrete 
systems subject only to controllability conditions on the resulting 
augmented system. The difference between the continuous and discrete 
augmented systems arises from the fact that the r repeated zero 
eigenvalues in the continuous-time domain’ map into r unity eigen- 


values in the discrete-time domain. 
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CHAPTER FOUR 


SIMULATION RESULTS FOR A DOUBLE EFFECT EVAPORATOR 


4,1 INTRODUCTION 
This chapter presents simulation studies of the eigenvalue 
assignment techniques presented in Chapters Two and Three applied to 
a model of the pilot plant double effect evaporator in this department. 
A 5-state linear discrete model is used in the simulation studies. 
The simulation results are divided into two main parts: proportional 
and PI control of the process, and are presented in this order. Before 
designing the controllers, a modal analysis of the evaporator model is 
performed to reveal the modal characteristics of the process model. 
The first part of the simulation study concerns the design 
of the proportional controllers and also discusses the performance 
of the resulting closed-loop systems. The recursive design technique 
described in Chapter Two is used and the open-loop eigenvalues are 
moved one by one. This control technique was selected for three 
reasons: 1) the evaporator model is "derogatory" as shown in 
Section 4.3, 2) eigenvectors at each recursive stage can be easily 
updated following the a eit outlined in Chapter Two, and 3) a 
better insight into the resulting process dynamics and controller 
gains can be attained by changing eigenvalues in a sequential manner. 
In the second part, PI controllers are designed such that 
eigenvalues of the original system are moved first, and then the 
repeated unity eigenvalues introduced by integral action are shifted 


using the methods presented in Chapter Three. 
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The performance of the controllers derived from the eigen- 
value assignment technique is compared with that obtained from optimal 
control theory to reveal the characteristics of the eigenvalue assign- 
ment techniques. Applications of optimal control theory to the 
evaporator were investigated extensively by Newell [22] using the 
five state linear model. The optimal feedback control laws used by 
Newell were derived using dynamic programming and minimizing a suitably 
chosen quadratic performance index. The discrete optimal feedback 
control laws used in this chapter were calculated in a similar fashion 
using the computer program, GEMSCOPE [47]. 

Since the eigenvalue assignment techniques dealt with in 
this thesis use the modal information of the system, it was necessary 
to accurately calculate the eigenvalues and eigenvectors of the open- 
loop system. At the present time the QR double iterative method seems 
to give the most accurate eigenvalues and eigenvectors of a general 
square matrix [45]. A program, CS201A [20], based on this QR double 


iterative method was used in this investigation. 


4.2 MATHEMATICAL MODEL 

A simplified schematic diagram of the pilot plant double 
effect evaporator is shown in Figure 4.1. The first effect is a 
natural circulation calandria type unit with a nominal feed rate of 
5 1lb/min of three percent triethylene glycol by weight. This feed 
is heated by a nominal 2 .1b/min of fresh steam. The second effect 
is an externally forced-circulation long-tube vertical unit with 
three 1" x 6' tubes and it utilized the vapor from the first effect 


as its heating medium. The vapor from the second effect is totally 
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condensed and the final product is about 1.5 1b/min of ten percent 
glycol. Tight vacuum control on the second effect maintains the 
necessary temperature differential across the tube wall. A more 
detailed description of the evaporator can be found in the thesis 
by Newell [22]. 

Models of the evaporator have been extensively investigated 
by Andre [2], Newell [22] and Wilson [48]. The model used in this 
thesis was derived by Wilson [48] by linearizing the nonlinear model. 
The process variables are employed in normalized perturbation form 
and the resulting model is a 5-state, linear, time-invariant model 


in the form of Equation (2-4). 
Z(t) = Ax(t)+Bu(t)+Dd(t) (2-4) 


Numerical values of matrices A, B and D can be found in Wilson [48] 
and the state vector x(t), control vector u(t), and disturbance 


vector d(t) are defined as: 


x, (t) W, (t) 
| x, (t) C(t) 
LE x(t) n hy Ct) 
x, (t) W, (t) 
x(t) C,(t) 
u, (t) SiC) 

BCE) oS au Ce) = Bi (t) (454) 
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where Wy> C > S, etc. denote normalized perturbation variable, 
(W Wiis ¢) 
as 1 
e.g. W os Tier (ss = steady state) 
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and the symbols are defined in the Nomenclature section. 


For the digital algorithms used in this study, the continuous 
model of the evaporator was discretized using the analytical solution 


to give the discrete model of the form of Equation (2-33), e.g., 
xICKEE) To OCT yack) oa ACE) uCkt) at OGL) dC kT) » «. (2-33) 


The numerical values of matrices $(T), A(T) and O6(T) are given 
in Table 4.1 for a sampling time T of 64 seconds. 

To design a PI controller for the evaporator using eigen- 
value assignment techniques, the eerie Z. in Equation (3-1) must 
be selected properly. At most 3 states can be selected for integral 
control due to the rank condition given in Section 3.3. Although 
arbitrary selection of the three variables is possible from among the 
five state variables, the holdups of the first and second effects 
Wy and Wo > 


on the physical importance of regulating these three variables. 


and the product concentration, Co> were chosen based 


Consequently, the numerical value of matrix i defined in Equation 
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Numerical Values of Matrices 


o(T) = 0.0 


-0.0119 
0.0116 


A(T) 


i} 
© 


. 1568 


-0.0137 


0.0137 


TABLE 4.1 


o(T) , 


T = 64 seconds 


0008 


9223 


0042 


- 0009 


0594 


-0.0817 


-0.0432 


On O92 


0.0871 


0.4377 


-0.1051 


0.1048 


-0.0406 


-0.0050 


0.0049 


0.0662 


-0.0058 


0.0058 


A(T) 


and 6(T) 


0.0 


0.0 
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(3-1) is given by 


(4-2) 
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The continuous-time augmented system can be formulated 
directly by combining Equations (2-4) and (4-1) and results in a eight- 
State variable model. As shown in Section 3.3.2, the discrete augmented 


model of the evaporator is expressed in the form of Equation (3-35) 
x[(k+1)T] = $(T)x(kT) + A(T) uCkT) + Q(T)d(kT) . (3-35) 


Although several approaches can be employed for the calculation of 
numerical values of matrices o(T), A(T) and eG), Newell and 
Fisher's approach was used in this study. Numerical values of 


matrices a(r).. A(T) and 6(T) are given in Table 4.2. 


4.3 MODAL CHARACTERISTICS OF THE EVAPORATOR MODEL 

The first step in the application of eigenvalue assignment 
techniques presented in Chapters Two and Three to the process control 
is to calculate eigenvalues and eigenvectors of the open-loop system 
model. For the discrete evaporator model, a modal analysis of the 


open-loop system gives: 
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Left Eigenvector Matrix: 


ON9871 0.0 0.0 0.6866 0.0 


-0.0014 0.9842 0.0087 -0.0006 0.6804 


V [vj >- Ve] =OR26020% O81771°" 0.999990 2393 BOR 2514 
0.0 0.0 0.0 0.6866 0.0 
0.0 0.0 0.0 0.0023 0.6884 
Right Eigenvector Matrix: 
1.013 0.0 ORU6254"0S0 0.0 
0.0 P-O1L7 68-0. T8024 "050 0.0 
W = [w,>- we] = 0.0 =020088* "00TS +40. 0 0.0 


—#5013 0.0012 0.1872 1.4565 -0.0048 


0.0 ~£50026 90018764 +070 £34553 


The above numbering of eigenvalues and eigenvectors is arbitrary, but 
convenient for the treatment that follows. However, vectors ia and 
w, inmatrices V and W are the left and right eigenvectors 

; ; siti ; 
associated with the i eigenvalue de> respectively. The elements 
of vector, vy, are normalized such that SV Suia I ue EV ier Gt i 
and <v,,w,> = 0 CO Ow 

It is clear that the evaporator model has two repeated 

eigenvalues of unity, which result in unstable open-loop responses 
for step changes in some of the inputs. Furthermore, the eigenvectors 
associated with these unity eigenvalues are linearly independent, thus 


the nature of the system is derogatory. It should be noted that eigen- 


vectors associated with the repeated eigenvalues are not unique, since 
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any vector resulting from a linear combination of these two eigenvectors 
also satisfies the definition of the eigenvector. However, the 
resulting two eigenvectors must be linearly independent. 

The diagonal elements of matrices Vand W are much 
larger in magnitude than the off-diagonal elements except for two 
columns, where one other element is comparable to the diagonal element. 
Since matrices o(T) and W are diagonally dominant, it is possible 
to associate each eigenvalue Ay with a single state variable of the 
system, i.e., the first eigenvalue AL is associated with the first 


state etc. The structure of the right eigenvectors reveals an 


X1> 
interesting fact that the three states corresponding to the stable 
eigenvalues of the system are completely decoupled from the two 
unstable eigenvalues, hence these three states will be stabilized 

in the presence of sustained disturbance. The two unstable eigenvalues 
interact in such a way that, in the presence of sustained disturbances, 
the first state Wy will exhibit an unbounded response independent 

of the fourth state W but the effect of a disturbance on the 


2? 
latter may be compensated partially by the interaction from the 
eigenvalue associated with the first state. This interaction of the 
two unstable eigenvalues on the fourth state will depend upon the 
type of disturbance. 

Figure 4.2 shows the simulated open-loop response of four 
state variables to a step disturbance of 20% in feed flow rate. The 
reasoning mentioned above explains the slow increase in the second 


effect level in contrast with the rapid, unstable response of the 


first effect level. The other two states, the first effect con- 
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centration and product concentration, eventually settle out to new 
steady state values, as expected. 

In order to stablilize the process, the two unity eigen- 
values must be moved to locations inside the unit circle in the 
complex plane. The mode controllability matrix H defined by 


Equation (2-10) was calculated to be, 


-0.0369 -0.0806 0.0 
0.0392 0.0 0.0 
H = 0.1569 0.0 0.0 (4-3) 
-0.0551 0.0020 -0.0279 


0.0567 ~1=0.0298., 10.0 


The ia row of matrix H represents the influence of the control 
vector upon the ayl eigenvalue of the system. However, since the 
system is derogatory with two repeated eigenvalues, the first and 
fourth rows of matrix H must be linearly independent for the system 
to be controllable [41]. The linear independence of the first and 
fourth rows is evident and the other rows have at least one non-zero 
element; thus the system is completely state controllable and all the 
eigenvalues can be shifted to arbitrary locations. It should be kept 
in mind that these two repeated eigenvalues cannot be changed by 

the simultaneous design mentioned in Section 2.3.2 and thus, at 


least, two stages of recursion are required [41]. 


4.4 PROPORTIONAL CONTROL OF THE EVAPORATOR 


The simulation results presented in this section include the 
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proportional controllers designed by assigning selected eigenvalues 

to the closed-loop system and also include the responses of the closed- 
loop system. Several responses are shown in the same figure in order 
to illustrate the effect of a particular design option upon the 
dynamics of the closed-loop system. Comparison of the performance of 
the controllers is mainly based on the magnitude of the gain elements 
and the dynamic responses of the closed-loop systems to typical step 
disturbances. In view of their physical importance in maintaining 
satisfactory control of the evaporator, the dynamic behavior of the 
first and second effect levels and the product concentration are 


emphasized. 


4.4.1 Separate Use of the Controls 

As mentioned in Section 2.3.3, the controls can be used 
separately to shift eigenyalues by specifying only one non-zero element 
in the vector g5 at each recursive stage. This results in the 
pairing of open-loop eigenvalues with the controls. 

For a multi-input system, the control matrix which assigns 
the desired eigenvalues to the closed-loop system is not unique and, 
in general, depends on: 

i) the pairing of open-loop eigenvalues with the controls, 

ii) the sequence in which open-loop eigenvalues are changed, 
and 

iii) the pairing of open-loop eigenvalues with the desired 
closed-loop eigenvalues. 
The effect of each of these design options is considered below. 


Since all five eigenvalues of 9(T) are controllable, 
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attempts were made to assign the set of desired eigenvalues, 


{p,} = {0.1, 0.2, 0.3, 0.4, 0.5}, to the closed-loop system. 


Effect of Pairing Open-Loop Eigenvalues with Controls 

In order to isolate the effect of pairing open-loop eigen- 
values with controls, the pairing of open-loop eigenvalues with desired 
closed-loop eigenvalues and the sequence in which open-loop eigen- 
values were shifted, were fixed. Thus, the Has element of the set 
{A 3° ye ie des ro} was moved to the Ae element of the set 
10i5:,80.2, 0.8 025, 0.4) at the i? recursive stage. 

In Table 4.3, the feedback matrices from Runs 1-3 show the 
effect of pairing eigenvalues with controls on the magnitude of the 
controller gain elements. The gain elements depend upon the pairing 
employed and change drastically from one run to another. In Figures 
4.3a and 4.3b, closed-loop responses to a step distrubance of 202 
increase in feed flow rate are shown using the three controllers from 
Runs 1-3. (The symbols used in the figure captions are defined in 
the Nomenclature section.) 

The pairing of eigenvalues with controls used in Run 1 was 
selected by inspecting the mode controllability matrix H and 
deciding which control has the largest influence on a particular 
eigenvalue. For example, from the numerical values of the elements 
in matrix H Ay is most influenced by Us» do by u)> and 
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dg by u,: 


with u, and Uy to avoid excessive use of ul: These pairings 


agreed with physical intuition except. for pairing ds with U,> 


The other two eigenvalues, hy and Aes were paired 


However, the resulting feedback matrix is not satisfactory since 
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some gain elements are excessively large and a highly undesirable 
transient results with large steady state errors (cf., Figure 4.3a). 
The pairing of open-loop eigenvalues with controls used in Runs 2 
and 3. was obtained by trial and error in order to ensure that the 
magnitude of the gain elements was reduced. A drastic reduction in 
the magnitude of gain elements was achieved in Run 3. The pairings 
used successfully in Run 3 seem to contradict the influence of controls 
on the open-loop eigenvalues shown by the numerical values of matrix 
H. For example, configuration used in Run 3 shows that ST which is 
associated with the first effect holdup was moved by manipulating 
the product flow rate. This phenomenon can be explained by the 
changes in the left eigenvectors associated with the unchanged eigen- 
values at each recursive step. Hence, the numerical value of matrix 
H also changes at each stage of recursion. In extreme cases, a 
controllable eigenvalue-control pair at the previous recursion could 
no longer be controllable at the next stage of recursion. 

Figures 4.3a and 4.3b show that if the closed-loop system 
has the same set of eigenvalues, then a better dynamic response to 
a feed flow disturbance is obtained as the magnitude of the gain 
elements decreases. Comparison of the feedback gain matrices for 
Runs 1-3 reveals the tendency that control effect is mainly achieved 
by.cu in Run 1 and u, in Run 2, while the control efforts are 


vs 


more uniformly distributed over all controls in Run 3. 


Effect of the Sequence of Changing Eigenvalues 


In order to study the effect of the sequence used in 


changing open-loop eigenvalues, the pairing of controls with open- 
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loop eigenvalues and the pairing of closed-loop and open-loop eigen- 
values were fixed. These pairings were the same as used in Run 3. 
The configurations of Runs 3-5 in Table 4.3 show the effect of the 
sequence of changing open-loop eigenvalues. The resulting feedback 
matrices are included to show their structure and the magnitude of 
the gain elements. The corresponding closed-loop responses to the 
20% step increase in feed flow rate are shown in Figure 4.4. 

The response of the closed-loop systems in Figure 4.4 shows 
that if the closed-loop eigenvalues are the same, better dynamic 
behavior will result for a controller with smaller gain elements. 
Steady state errors in the state variables due to three different 
disturbances are shown in Table 4.4 for Runs 3-5. Although, the 
magnitude of steady state errors varies remarkably from one distur- 
bance to another, the same trend, that the smaller the gain magnitude, 
the better the dynamic behavior, seems to hold for other disturbances. 
Furthermore, it is observed from Table 4.4 that the feed concentration 
disturbance has the most significant effect in comparison with the 
other disturbances, while the same disturbance has the least effect 
on the open-loop resonse. 

The feedback matrices for Runs 3-5 show that the matnitude 
of gain elements also depends on the sequence in which open-loop 
eigenvalues are changed. It is observed from Runs 1-5 that the control 
effort should be distributed over the available controls to obtain 
lower gains and good dynamic behavior. It should also be noted that 
a particular eigenvalue-control pair in one sequence may not be 


possible in a different sequence of eigenvalue shifts due to a change 


in left eigenvector directions. 
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TABLE 4.4 


Steady State Errors for the Closed-Loop Systems - I 


; Steady State Beroce< to 
Run No. State Variables NB epi ee ee 
T2040 = 207 (CP 2 e204 eTF 
Wi 1.28 0.01 =a 8 ses 
C, U.Ol =L9 97 0.0 
3 hy 5.68 Org 0.0 
W. Paps S88 Geld =0630 
C., DA 8 rae BS JS 0.01 
Wy Vege 3.36 = 0.32 
C, 05605) 7=19' 5/2 020 
4 hy Die Ol, Ons L 0.0 
W, She Os) Bae - 0.36 
C, JAA ARCS RS es Oe) 0.01 
W, = 6 539 85.80 ork 
Cc) =—10'.98 AY CA 25 
5 hy 4.72 9.08 Ould 
Wo O25 84.84 =e MATT 
C, = 5.86 45.35 de dal 
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* ° 
Steady State Errors were taken at t = 60 min and are expressed as 


a percent of the initial states. 
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Effect of Pairing Open-Loop and Closed-Loop Eigenvalues 


The configurations in Runs 3, 6 and 7 are different only 
in the pairings of open-loop and desired closed-loop eigenvalues that 
were used. Although the magnitude of gain elements changes with 
variations in the open-loop and closed-loop eigenvalue pairings, this 
effect is small compared with the other factors previously discussed 
and the resulting feedback matrices have the same general features. 
The effects of pairing open-loop and closed-loop eigenvalues on the 
dynamic response are seen in Figure 4.5 for a 20% step disturbance 
in feed flow rate. The difference in dynamic behavior is very small. 
In conclusion, it is observed that the dynamic behavior 
of a closed-loop system depends highly on the configuration used in 
shifting open-loop eigenvalues. Of the three options discussed, the 
sequence of changing eigenvalues and the pairing of open-loop eigen- 
values with controls have significant effects on both the gain magni- 
tudes and the closed-loop dynamics. If the closed-loop eigenvalues are 
the same, the closed-loop response becomes more satisfactory as the 
gain elements decrease in magnitude. It is interesting that, in 
the resulting feedback control systems, the holdups of the first 
and second effects govern the steam (u,) and product (ug) flow 
rates while the bottom flow rate (u,) from the first effect is 
governed mainly by the first and second effect concentrations. Further- 
more, the controls, in general, are not significantly affected by 
the enthalpy of the Firat effect, as indicated by the relatively 


small elements in the third column of KEB in Table 4.3. 
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4.4.2 Simultaneous Use of Controls 

If all the elements in vector 85 in Equation (2-21) are 
specified to be non-zero, then all controls are used to change an 
eigenvalue at each recursive Stage. An arbitrary selection of 85° 
which gives non-zero elements and preserves the controllability of 
the eigenvalue to be changed, would result in a controller matrix 
which assigns a desired set of eigenvalues to the closed-loop system. 
Although Equation (2-32) is valid for gain minimization when a single 
eigenvalue is changed, vector 8, chosen by Equation (2-32) at each 
recursive stage will guarantee the controllability of an eigenvalue 
to be shifted but the resulting controller gains will not necessarily 
be minimal. 

In this section, open-loop eigenvalues are changed recur- 
sively with 85 specified at each stage by Equation (2-32). The 
effects of 

i) the sequence of changing open-loop Bicenvaticess and 

ii) the pairing of open-loop and closed-loop eigenvalues 
on the feedback matrix and the dynamic response are investigated. 
Furthermore, a controller designed using an arbitrary set of vectors 


{g.}, is compared with one of the controllers designed using 


Equation (2-32). 
The set of desired closed-loop eigenvalues chosen in this 
section is the same as was used in previous subsections. The 


configuration used in each design and the resulting feedback matrices 


are summarized in Table 4.5. 
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Effect of the Sequence of Changing Open-loop Eigenvalues 

In Runs 8-10 in Table 4.5 and Figure 4.6, the pairing of 
open-loop eigenvalues with desired eigenvalues was fixed and the 
sequence in which open-loop eigenvalues were changed, was varied. 
As in the case of using controls separately to shift eigenvalues, 
this effect is significant for both the controller gains and the 
dynamic behavior of the resulting closed-loop system. It can readily 
be seen from Figure 4.6 that the controller resulted from Run 10 
shows the best response with the smallest offsets in the four state 
variables. In Table 4.6, steady state errors due to typical step 
disturbances are given for the closed-loop systems resulted from 
Runs 8-10. The same dependence of the closed-loop system performance 
on the type of disturbance as was observed in the separate use of 
controls can be seen from Table 4.6, Although individual states 
exhibit different steady state errors depending on the type of 
disturbance, in general, the dynamic behavior of OU aigcmdeiogs 
system is better when the magnitude of the gain elements is smaller 


for the same closed-loop eigenvalues. 


Effect of Pairing Open-Loop and Closed-Loop Eigenvalues 


The feedback matrices obtained from Runs 10-12 in Table 4.5 
show that the pairing of open-loop and closed-loop eigenvalues can 
effect the magnitude of the gain elements. Figure 4.7 compares the 
dynamic response of these three closed-loop systems to the 204 increase 
in feed flow rate. In contrast with the case of using controls 
separately, the gain elements change significantly depending on how 


the open-loop and closed-loop eigenvalues are paired. However, the 
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Figure 4.6 Simulated Effect of the Sequence of Changin 
Use of Controls (+20% F/P) 
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TABLE 4.6 


Steady State Errors for the Closed-Loop Systems - II © 


Steady State Errors to 


State Variables 


+20Z% F -20% CF +20% TF 
W, ~1.80 36.38 0.38 
C, 02108) 3-15) 62 0.09 
h, Sey: Dace 0.05 
W, Male h Gece WAS OW amet a’ 
es, Hi 09m = 0.72 0.39 
W 22209 45.81 0.58 
ae =O8 08 4y 1.119530 0.14 
hy 5.26 AVON 0.08 
W, Bee T ay aang ee | Nado7 
oy #1 89 9.57 0.60 
Ws 0.85 (2-06) v0.23 
oH DS ae yt 0.01 
hy 5.58 0.95 0.01 
W, 1.00 P107en 0-14 
C 0.46" 14.36 0.05 
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controller matrices are similar and the dynamic behavior of the three 
closed-loop systems is almost indistinguishable. The effect of pair- 
ing open-loop and closed-loop eigenvalues is far less significant than 


the effect of changing the order in which open-loop eigenvalue are shifted. 


Effect of the Choice of a Vector Set, tg 

In Runs 13 of Table 4.5, the design configuration is identical 
to that of Run 10 except that the vector set, {g,}, was specified 
arbitrarily instead of using Equation (2-32). The resulting set of 


vectors, {g,J, are as follows, 


ui an 2 
g, = i 5 &> = = ih . 83 = il ; 
it it 1 
iL it 
Ws 2 ‘ g5 = ib 
uh Z 


Comparison of the controller matrices from Runs 10 and 13 and the 
transient peaaonee in Figure 4.8 indicates that excessively high gains 
in the feedback controller can be avoided by making use of Equation 
(2-32) in the selection of 85° 

The results of this section can be summarized as follows: 

i) As in the case of using controls separately, better 
dynamic response is expected as the magnitude of gain elements is 
reduced if the closed-loop system have the same eigenvalues. 


ii) The sequence of changing eigenvalues has a more signifi- 
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cant effect than the pairing of open-loop and closed-loop eigenvalues 
on both the controller gains and the closed-loop dynamics. 

iii) The pairing of open-loop and closed-loop eigenvalues 
has a significant influence on the controller gains, but the dynamic 
responses of the resulting closed-loop systems are only slightly 
affected. 

iv) Equation (2-32) can be utilized in specifying a set 
of vector, {g,} in order to avoid extremely large controller gains, 
which might occur when a set of arbitrary vector {g.} are used. 
Although, each vector 8, selected by Equation (2-32) locally 
minimizes the gain matrix at the hee stage, there is no guarantee 


that the final gain matrix is minimized. 


4.4.3 Effect of the Desired Closed-Loop Eigenvalues 

Another important question, independent of the design options 
discussed in the previous two subsections, is how the magnitude of 
the closed-loop eigenvalues affects the dynamics of the closed-loop 
system and the feedback gain matrix. As an attempt to investigate 
this design choice, two different sets of desired eigenvalues, namely 
{0':2, O44, he @unchanged)s 0.5, O.6} wand) 40:4,6 0.54 0..6,.0.7, 0.5835 
were assigned to the closed-loop systems in Runs 14 and 15 respectively. 
The design configurations used in Runs 14 and 15 and the resulting 
gain matrices are given in Table 4.7 together with those of Run 10 
for comparison. The simulated responses of the resulting closed- 
loop systems are shown in Figures 4.9-4.11 for step disturbances in 


feed flow rate (+20%), feed concentration (-20%) and feed temperature 
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The feedback matrices in Table 4.7 show that the largest 
gain element Scere sees significantly as the magnitude of the closed- 
loop eligedvataae is increased. This is an expected result since less 
control will be required as the open-loop eigenvalues are moved smaller 
distances. The dynamic behavior of the closed-loop system to the dis- 
turbance in feed flow rate, as shown in Figure 4.9, is not surprising 
since the closed-loop eigenvalues and offsets are smallest in Run 10. 
However, in Figure 4.10 where a feed concentration disturbance is 
considered, the behavior observed in Figure 4.9 is completely reversed. 
This behavior is hardly expected judging solely from the mangitude of 
the closed-loop eigenvalues. The effect of a feed temperature dis- 
turbance on the closed-loop dynamics is least significant, as shown 
in Figure 4.11, compared with the other two disturbances. An interest- 
ing observation is that Runs 14 and 10 yield almost identical responses 
to the step disturbances. Comparing these dynamic responses and the 
magnitude of the gain elements, the controller from Run 14 performs 
better than that from Run 10. However, the new steady state is 
reached faster as the magnitude of the closed-loop eigenvalues is 
reduced, as would be expected. 

From these abienes tionee it is clear that a set of smaller 
closed-loop eigenvalues does not necessarily yield better performance 
of the closed-loop system. It is observed that the dynamic response 
is highly disturbance dependent as was the case in the previous two 
subsections. This behavior can be explained qualitatively in terms 
of the directions of the resulting closed-loop eigenvectors. If 


the resulting closed-loop left eigenvectors happen to be oriented in 


‘end os linia M AP aA 
abseats. on3 ¥ ob gaat i an 
esol Sonbe tluesy Roaseas wae et | 8 

sehbens bayom oré -esolebacgie goot-noge aa an 4 
SED) asta 03 mereye qoot=badets oid: to solaiied : 
gut cdtirs ton ek .e@.e baeaiis) pi owede ,9ae3 woth sd 
; OL twa ot sesileme SB etree seulement sere i 


 beBtovet Leia st R.s sagged of bowabeds 1 


to sbusignsm sat sox? atiegpi aniapet basoeane. ‘eibeed 8 


~elb siutsvegsss bee% 6 ‘Fo! dogtts: SAP te. 

= words a6 japoriingts se801 at eo taeayh qoot-be : 9 etd 96 
~Seeretal oh _seonkdsuseth: ows — sits dat borage abe: 7 ones 
ssumoqes igotsnsbh- teomis phage OL bine AL enue dei ‘ak -nolsnrsegdo ga be yy 
add bas euenogas otmanyb seeds antregaod -eeoandsus0t qt os 03 


amiotrey al out mor isi Loxtao9 sig “ 


et odeda eres wan alt cadieall a j oa sot ‘nats nods 
et, 2oulayasgis oot-baaols oul bide sig ht sek 


wollen to 192 § iad seals, at ok pao seus 
spadaretrag reused bioby eilxsessoon 308 aod 
enniodaey pie. ef ted Devrronde wb: oi 


: wa ayotvers re HE ae59 ‘its aw eB ‘ 


eats bia xfoviseatieap bontalgqes ad a9 


ce eroaavangts rein * 


= 


87 


the same direction, then from Equation (2-7), it can be shown that 
the components of the right eigenvectors are very large. Hence, 
those disturbances whose corresponding coefficient vectors (columns 
of matrix OCT) are relatively colinear to the left eigenvectors 
will have significant effects on the dynamic behavior of the state 
variables, even for small eigenvalues. This point will become more 
apparent in the next subsection where an optimal control law is 


compared with control laws derived eigenvalue assignment techniques. 


4.4.4 Comparison With an Optimal Controller 

Since controllers derived from optimal control theory are 
not explicitly concerned with shifting eigenvalues and the controllers 
derived from eigenvalue assignment techniques are not unique for a 
multi-input system, a meaningful comparison of these two types of 
controllers is difficult. However, if the closed-loop systems 
derived from optimal control theory and eigenvalue assignment 
techniques have the same set of eigenvalues, a comparison of closed- 
loop dynamics will reveal the effect of the eigenvector orientations 
of the closed-loop systems. 

The Suititnall control law used in this study was calculated 
from the GEMSCOPE computer program [47], and is identical to one of 
the controllers used by Wilson [48]. In view of the physical 
importance of the state variables to be controlled, the state weight- 
ing matrix, Q = diag [10, 1, 1, 10, 100], used in previous studies 
[22, 48] was chosen. The control weighting matrix, R, was set 
equal to the zero matrix. The eigenvalues of the closed-loop system 


using this optimal control law were calculated to be {0.9002, 0.2706, 
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=e. woe. 1x on 4 x Von ay This spectrum of eigenvalues is quite 
interesting since three eigenvalues are extremely small, one is 
relatively small and the other one is large in magnitude. This set. 
of closed-loop eigenvalues was used as the desired set of eigenvalues 
in designing controllers by eigenvalue assignment techniques. 

In Table 4.8, feedback gain matrices derived from optimal 
control theory and two eigenvalue assignment techniques are shown 
along with the closed-loop eigenvalues and the design configurations 
that were employed. The controllers derived in Runs 16 and 17 are 
those with the smallest gain magnitudes that were obtained by trial 
and error; the largest gain elements are significantly smaller than 
the comparable elements in the optimal controller. A pair of complex 
eigenvalues was obtained in both Run 16 and Run 17. This creation of 
a pair of complex eigenvalues always occured when more than two open- 
loop eigenvalues were to be moved to very small values (~0). During 
this study, it was observed that the closed-loop eigenvalues were 
very sensitive to the numerical accuracy of the left eigenvector used 
at each recursive stage. This is especially true when some of the 
desired eigenvalues are near the origin and the word length of the 
digital computer is limited. For this reason, the creation of the 
complex eigenvalue pair could not be eliminated even by using "extended 
precision" on the IBM 360 computer. 

Figures (4.12)-(4.14) compare the closed-loop responses for 
the optimal controller and controllers derived using two eigenvalue 


assignment techniques for step disturbances in feed flow rate, feed 


concentration and feed temperature. The optimal controller effectively 
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reduces the influence of three different types of disturbances on 

the heavily weighted states while the performance of the other two 
controllers is still dependent on the type of disturbance present. 
Since the closed-loop eigenvalues are almost identical for these 

three closed-loop systems, the different dynamic responses are entirely 
due to the orientation of the closed-loop eigenvectors. However, 
direct comparison of the closed-loop eigenstructure is impossible 
because of the pair of complex closed-loop eigenvalues resulting from 
Runs 16 and 17. Instead, the eigenstructure of the closed-loop 

system obtained from Run 10 is compared with that of poe controller 
to see the effect of the closed-loop eigenvector orientation. Although, 
the resulting closed-loop eigenvalues from Run 10 and the optimal 
controller are quite different individually, the magnitudes of gain 
elements are comparable and the closed-loop dynamics obtained from 
Runs 10, 16 and 17 show the same disturbance dependence. Furthermore, 
it will be seen that the size of the closed-loop eigenvalues is less 
important than the orientation of the closed-loop eigenvectors in 
determining the dynamic behavior of the closed-loop system. 

In Tables 4.9 and 4.10, the eigenstructure of the closed- 
loop systems resulting from Run 10 and the optimal controller are 
presented along with the mode disturbance matrices. The eigenvectors 
in Tables 4.9 and 4.10 are normalized such that Equation (2-7) holds 
and the Euclidean norm of each left eigenvector is unity. Then the 
determinant of the left eigenvector matrix is an indication of how 
well distributed the eigenvectors of the closed-loop systems are in 
the state space. The effect of poorly distributed left eigenvectors 


is directly reflected in the magnitude of the elements in the right 
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Eigenstructure and Mode Disturbance Matrix of the 
Closed-Loop System of Run 10 
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Table 4.9 - continued 
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Eigenstructure and Mode Disturbance Matrix of the Closed-Loop 
System from the Optimal Controller 
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Table 4.10 - continued 


Mode Disturbance Matrix, 9 
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eigenvector matrix. The comparison of the numerical values of these 
determinants of the closed-loop systems derived from Run 10 and the 
optimal controller with the value of 0.45 for the open-loop system 
confirms the magnitude of elements shown in the respective right 
eigenvector matrix. 

Eivar zero initial) statejwe. 3, x(0) t= Q and a constant 
disturbance vector, e.g., d(k) = d is assumed, the state at the 
ee sampling time is related to the eigenstructure of the closed- 


loop system (with distinct eigenvalues) in the following manner: 


x(kT) = WACKY" 9d 
or 
ae wat 
x(kT) = ) <8 v,,d>p, (kw, (4-4) 
i=l 
where 
A(k) = diag [o, (kK) 0 (k),--- 50 Ck) ] 
and 
1-p 
= Ail wea 
p, (k) 13, ( n) 
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which has small components along some of the column vectors of the 
disturbance matrix 8, then the effect of corresponding disturbances 

on the states will be insignificant. But, for other disturbances, the 
response of the closed-loop system will suffer from the large elements 
in the right eigenvector matrix. For the closed-loop system of Run 10, 
Table 4.9 shows that the first four left eigenvectors are clustered 
around a direction which has a significant component along the direction 
of the second column vector of matrix 9. Only the last left eigen- 
vector is fairly well separated from the other left eigenvectors. Con- 
sequently, the elements of the last right eigenvector are very small 
compared with those of the other right eigenvectors. This explains 

why the closed-loop response of Run 10 is poor for the feed con- 
centration disturbance despite the small eigenvalues. 

On the other hand, the left eigenvectors of the closed-loop 
system resulting from the optimal control law are fairly well dis- 
tributed and prevent excessively large elements in te right eigen- 
vector matrix. Considering the comparable magnitude of the elements 
in matrix ov for the cases of Run 10 and optimal controller, the 
difference in the closed-loop dynamics is mainly to the directional 
effect of the right eigenvectors in the state space. Furthermore, 
although the optimal controller was obtained without explicit reference 
to the closed-loop system, the resulting closed-loop eigenstructure 
is such that the effects of the disturbances on the heavily weighted 
states are minimal. This is achieved not only be the small closed- 
loop eigenvalues but by eigenvector directions such that the combined 


effect of eigenvalues and the associated eigenvectors rejects the 


influence of the disturbances or non-zero initial state as much as 
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possible. However, even the optimal controller seems to be unable 

to regulate the heavily weighted states to the same extent depending 

on the disturbance. In Figure 4.12, despite the same weighting on 

the first and second effect holdups, it is seen that the response 

of the first effect level has a much larger offset than that of the 
second effect level. On the other hand, the first effect concentration, 
which is ten times less weighted than the level, exhibits smaller 
offset than the level. This leads to the conclusion that the direction 
of the closed-loop eigenvectors cannot be controlled arbitrarily, 

even if the closed-loop eigenvalues are not prespecified, when the 
system is to be stabilized and the number of controls is less than 

the number of states. 

The comparison of the optimally controlled system with the 
closed-loop systems derived by the eigenvalue assignment techniques 
reveals the importance of eigenvector directions of the closed-loop 
system in the regulation of the states. It was shown that a closed- 
loop system with well distributed eigenvectors in the state space 
exhibits better dynamic performance than a closed-loop system with 
poorly distributed eigenvectors. Unless the number of the controls 
is equal to the number of states, there seems to be a limitation in 


controlling closed-loop eigenvectors when the system is to be 


stabilized. 


4.5 PROPORTIONAL PLUS INTEGRAL CONTROL OF THE EVAPORATOR 
The simulation studies of the evaporator presented in 
Section 4.4 confirmed the well-known fact that, in the presence of 
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completely in any of the states using only a proportional controller. 
Hence, it is the purpose of this section to eliminate offsets in some 
of the states by using a PI controller designed by the eigenvalue 
assignment techniques presented in Chapter Three. 

The state variables chosen for the integral control action 
and the resulting augmented open-loop system were presented in 
Section 4.2. This section is concerned with the egrets of the 
following design factors: 

i) the effects of the two algorithms used in the design 
of integral portion of the PI controllers, 

ii) the effects of the design configurations used in the 
first step of the PI controller design, and 

iii) the effects of the different sets of the closed-loop 
eigenvalues assigned in the first step of the PI controller design, 
on the performance of the resulting systems. 

As a reference controller for comparing ae perience 
of the resulting closed-loop systems, an optimal PI controller was 
designed using the augmented open-loop system. Due to the increased 
dimension of the augmented system, a state weighting matrix, Q, of 
dimension (nt+tr) was apeeaeied in the performance index. Further- 


more, to reduce the magnitude of the gain elements, a non-zero control, 


weighting matrix, R, was selected. The numerical values of ‘the 
matrices, 6) and R, used in the design of an optimal controller 
are given by; 
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The resulting optimal PI controller is given in Table*4.12. i. The 


corresponding closed-loop eigenvalues are, 
{p,} =-10'.1299° 40 .sap6+0612491, 036566, 0.97549, 0.7824, 0.8980, 0.9000} . 


In the design of the PI controllers using eigenvalue assign- 
ment techniques, these eigenvalues were specified as the desired closed-~ 
loop eigenvalues. However, since the optimal PI controller produced 
a pair of complex conjugate eigenvalues, the desired eigenvalue cor- 
responding to these complex eigenvalues were replaced by two real 
eigenvalues. This was done to avoid the difficulty of handling 
complex eigenvectors. 

The design configurations used in the eigenvalue assignment 
techniques are summarized in Table 4.11. The corresponding PI 


controllers are tabulated in Table 4.12. 


4.5.1 Comparison of the Recursive and Simultaneous Designs 

The design configurations used in Runs 18 and 19 differ 
only in the second step, where the three unity eigenvalues are shifted 
to the desired locations using recursive and simultaneous approaches, 
respectively. Controls were used simultaneously with the g. vectors 
calculated from Equation (2-32) at each stage of the first design 
step. The design configurations of Runs 20 and 21 are analogous to 
those of Runs 18 and 19, except that, in the first step, controls 
were used separately to shift the five eigenvalues of matrix o(T). 
The pairing of controls and eigenvalues that were used in Runs 20 and 


21 are shown in Table 4.11. 


Although the same design configurations were used in the 
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first steps of Runs 18 and 19, significantly different proportional 
feedback gains were obtained. Similarily, the proportional feedback 
matrices resulting from Runs 20 and 21 are quite different. However, 
the general features of the integral feedback matrix, including the 
Signs of the gain elements, seem to be preserved if the same design 
policy is used in the second step of the design. This can be seen 


from a comparison of the K. matrices resulting from Runs 18 and 20, 


L 
and Runs 19 and 21. 

The transient responses of the closed-loop systems using 
the PI controllers of Runs 18 and 19 are shown in Figures 4.15 to 
4.1/7 for step disturbances of +204 in F, CF and TF, respectively. 
The transient response for the optimal PI controller was also plotted 
for purposes of comparison. As in the case of proportional control, 

a disturbance in feed concentration has the most significant effect 
and a disturbance in feed temperature has only a negligible effect 

on the dynamic behavior of the closed-loop system. In Figures 4.18 
and 4.19, the closed-loop response using the PI controllers derived 
in Runs 20 and 21 were plotted for step disturbances in F and CF 
only, since the effect of the disturbance in TF was also negligible 
in these two runs. 

The transient response varies depending on not only the 
disturbance but also on the design method used in the second step 
of PI controller design. The comparison of transient responses in 
Figures 4.15 and 4.19 indicates that the new steady states are 
reached faster ean the simultaneous design is used in the second 


step of PL controller design than when the recursive design is used. 
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This trend is apparent in all the responses regardless of which 
disturbance is present. For both the design methods, a feed flow 
disturbance produces the same type of transient response, but the 
response to the feed concentration disturbance changed significantly 
depending on the design method used in the second step. However, 

for the product concentration, Co» which is the most important 
variable in this particular application, the Sea ea ae design 
approach provides better transient responses (i.e., smaller dewiations 
from steady state) than the recursive approach (cf., Figures 4.16 


endi4019) . 


4.5.2 Effect of the Design Policy Used in the First Step 

The design configurations of Runs 18 and 20 differ only in 
the first design step and similarily for Runs 19 and 21. Hence, the 
difference between the closed-loop systems reflects the effect of 
the design approach used in the first step of the design. However, 
it is difficult to generalize about the effects on the feedback 
control matrices and the closed-loop transient responses. 

From Equation (3.56), the final expression for the propor- 
tional feedback matrix is given as, Kes = [K-TK 7 {1 +(gtQK-L) "1. 
Since this expression shows the explicit relationship between the K 
and K controller matrices designed in the first and second step, the 


it 


structure of KERB and the transient response of the closed-loop system 


depend significantly on the magnitude of the gain elements of matrix Kr 


4.5.3 Effect of the Desired Closed-Loop Eigenvalues 


In Run 22, only four eigenvalues of matrix 9(T) were 
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shifted to the locations given in Table 4.10. The other design options 
are the same as in Run 19. 

A comparison of the resulting feedback matrices for Runs 19 
and 22 does not reveal any common features between controllers. The 
Structure of the controller matrices is quite different when a dif- 
ferent set of eigenvalues are specified in the first step, even if 
the design configurations are identical in the second step. 

The transient response of the closed-loop system of Run 22 
is quite disturbance-dependent as was the case for other closed-loop 
systems. A disturbance in F results in a response similar to the 
one in Figure 4.15 and the response to the TF disturbance is very 
similar to Figure 4.17. For this reason, only the response to the 
CF disturbance is shown in Figure 4.20. It is apparent from 
Figure 4.20 that.the closed-loop eigenvalues assigned in the first 
design step have a significant effect on the shape of the transient 


response of the resuiting closed-loop system. 


4.5.4 Summary and Interpretation 

This summary is based on the performance of the closed-loop 
systems which were presented in this section. 

1) As expected from theoretical considerations, offsets 
were eliminated in three states, namely Wy> W. and C,; for all 
three constant disturbances. The elimination of offsets can also 
be shown analytically using the discrete representation of the 
augmented system but a rather lengthy derivation is required. 

2) For a given disturbance, each of the five state variables 


converged to the same steady-state value regardless of.which-Pi con- 
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troller was used. This phenomenon can be explained qualitatively 

in terms of degrees of freedom or in terms of the controllability 

of the augmented system. For example, the controllability condition 
indicates that only r (< m) state variables can be used in integral 
control. This means that only r_ state variables can be independently 
driven to arbitrary values and the remaining (n-r) state variables 
are then fixed since there are no more degrees of freedom available. 

3) The transient responses of the closed-loop systems 
resulting from PI controllers are highly dependent upon the distur- 
bance present. As observed in the study ae eh controls, 
the CF disturbance has the most significant effect on the closed- 
loop dynamics and the TF disturbance has a negligible effect. Since 
‘the recursive design method used in the first design step was identical 
to that used in the proportional controller design, this sensitivity 
to disturbances seems to be due to the recursive control policy used 
in the first step. 

4) The simultaneous design method tends a give better 
transient responses, namely, smaller deviations from the desired 
steady states and a faster recovery from disturbances. [In this 
study, a diagonal matrix consisting of the desired eigenvalues was 
used as matrix D in Equation (3-55). Clearly, the eigenvectors 
of matrix D are less colinear than those of the closed-loop 
rx vr matrix designed by the recursive design. Thus, it is expected 
that the resulting eigenvectors of the augmented system are also 
well-distributed when the simultaneous design is used. In view of 
the discussion given in Section 4.4, this could be a possible 


explanation for the better transient responses. 
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CHAPTER FIVE 


CONCLUSIONS 


The purpose of this work was to investigate the application 
of existing eigenvalue assignment techniques based on state feedback 
to the design of multivariable control systems. In particular, an 
evaluation of the design options available when the eigenvalue assign- 
ment was based on a modal analysis, was one of the main objectives 
of this study. The other objective was to extend eigenvalue assign- 
ment techniques to the design of multivariable PI controllers. 

A new theoretical result for PI control systems is that 
the integral control matrix must have full rank in order to arbi- 
trarily assign all of the eigenvalues of the closed-loop system. 

The constructive proof of the proposition provided the basis for two 
modified alogrithms for PI controller design by eigenvalue assign- 
ment techniques. | 

Proportional controllers were designed for a double effect 
evaporator using several recursive approaches and the resulting closed- 
loop dynamics were simulated. In general, the simulation results 
revealed that the performance of the closed-loop evaporator system 
was highly disturbance-dependent and not always satisfactory when 
these proportional controllers were used. The unsatisfactory behavior 
of the closed-loop system was mainly due to the orientation of closed- 
loop eigenvectors in the state space. It was shown in simultation 
studies that the PI controllers designed by eigenvalue assignment 


techniques were able to eliminate offsets in the state variables 
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Subjected to integral control. 


5.1 SUMMARY 

Since a rather detailed discussion of the simulation results 
for the evaporator system was given in each section of Chapter Four, 
only the general conclusions from the simulation studies are summarized 


here. 


Conclusions from the Proportional Control Studies 

The design options available in the recursive design methods 
have significant effects both on the dynamics of the closed-loop 
systems and the resulting controller matrices. In this regard, the 
sequence in which the open-loop eigenvalues are changed to the 
desired closed-loop eigenvalues had the most significant influence 
on the performance of the closed-loop system. However, effective 
utilization of these and the other design options was not obvious 
a priori and was generally achieved by trial and error. 

For the same design configuration and desired closed-loop 
eigenvalues, the closed-loop dynamics became more satisfactory as the 
largest feedback gain element was reduced. 

The new steady meee was reached faster when the closed- 
loop eigenvalues were small as would be expected, but the performance 
of the closed-loop system did not necessarily improve as the magnitude 
of the closed-loop eigenvalues became smaller. Furthermore, changing 
all the controllable eae eigenvalues to smaller values did not 
necessarily result in as improved dynamic response of the closed- 


loop system in comparison with the case where only some of the con- 


3 8 a Mi 
i 
ae 
i] ie f ' 
: i a, i - A 
q ae a ree} ' a fa uli f har 
: - } 


“ 
int Bk 
pa 


eiivess sotzalumia nid, bai noiseusath bottaiob ere x o on 
-auekt yesgeds io wigs’ wind ‘ak novig eae sow sot 33 2s 


jsstyanigne Sag a Side nob sutume aad 3 i 


pe ons a 
Betbure x orineD Lanokr to 


ehodtom ngbasb sy baraex atti a seating crops mateo pet. pes ive mM: 
qool-hagefo sity Be a hme oda ae, fied. Bide if: bag tog ; 

efit , brseox xT, -eoninem aloxano9, a Le pe Ke 
‘silt 09 bagenis on samdovnagie Solange women 
soueul tat gnanthing te, daom sii bad. Palen - ts 
ovdtostie .7evswoH asdiege qosi-bseots sacha we a 
enoivds tox egw 2eotsgo., agieab ai eal dl bail 
“tongs bins tabi’ oe | 

qool=bhaepis boriesb bas id | 


ot 2B 3 nodonteksee arom ages: 


~bexo arts farts sede stows aow, soe, : 
salinity al te avo 3 ; } 


sis it esctow + ssttaa OF) sevtevangie : } | 


120 


trollable open-loop eigenvalues were changed. Since the closed-loop 

dynamic behavior is determined by both eigenvalues and eigenvectors, 
manipulation of only eigenvalues is not always sufficient to improve 
the dynamic response of the closed-loop system. 

The closed-loop dynamics were very sensitive to the type 
of disturbance present. The feed concentration disturbance had the 
most significant effect on the dynamics of all the closed-loop systems 
while a feed temperature disturbance had very little effect. The 
disturbance-dependence of the closed-loop systems was due to the 
orientation of the closed-loop eigenvectors in the state space. 
Furthermore, the closed-loop eigenvectors tended to be oriented in 
the same direction in the state space. This could not be directly 
controlled by the designer using eigenvalue assignment techniques. 

In the evaporator simulation study, the closed-loop eigen- 
values assigned at each recursive step were very sensitive to the 
numerical accuracy of the updated eigenvectors when the magnitude 
of the desired eigenvalues was small. Hence, in some applications, 
the straightforward algorithm for updating eigenvectors may not be 
accurate enough to assign the desired eigenvalues, unless the 
numerical techniques are carefully evaluated. 

In conclusion, the simulation results revealed that in 
process control assigning eigenvalues of the closed-loop system is 
not sufficient and the closed-loop eigenvectors are also important 
to the same extent in ensuring satisfactory performance of the 
closed-loop system. The eigenvalue assignment techniques guarantee 
only the shift of open-loop eigenvalues to the desired locations but 


does not assign the closed-loop eigenvectors unless the available 
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design freedom can be exploited in some systematic way. However, 
the utilization of this design freedom will, in general, result in 
a non-linear optimization problem. Thus simplicity, which is the 


main advantage of eigenvalue assignment techniques, will be lost. 


Conclusions from the Proportional Plus Integral Control Studies 


The PI controllers designed by eigenvalue assignment 
techniques eliminated offsets in the properly chosen state variables 
of the evaporator system, regardless of the type of sustained dis- 
turbance. However, comparison of an optimal PI controller with PI 
controllers designed by eigenvalue assignment techniques showed that, 
although the same steady-state was reached, the transient response 
of the closed-loop system was more satisfactory when the optimal PI 
controller was used. Furthermore, for the PI controller derived from 
eigenvalue assignment techniques, the transient responses of the 
closed-loop systems depended on the types of disturbances present 
and the design configurations used. In view of the important role 
of the eigenvectors in the dynamic behavior of a system, this 
variation in the closed-loop dynamics was due to the eigenvector 
orientation, which depended on the design configurations used. 

For the two algorithms presented in Chapter Three, all 
the closed-loop systems showed the same disturbance-dependence. Since 
this reflects the initial design policy used in the first step, a 
satisfactory proportional controller, K, is more important than 
the integral control matrix, Ke if the two step design scheme is 
to be used. However, the simultaneous design approach used in the 


second step, in general, gave better transient responses of those 
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selected state variables than did the recursive approach. This again 
reflects the importance of the closed-loop eigenvectors in the 
dynamic behavior of the closed-loop system. 

In conclusion, the simplicity of eigenvalue assignment 
techniques can be advantageously used in eliminating offsets in 
selected states. For the two step design schemes used in this study, 
more emphasis should be directed to the design of the initial pro- 
portional control matrix, K, to get satisfactory transient 
responses. However, if the standard eigenvalue assignment techniques 
are applied directly to the augmented system, the relation between 
eigenvectors and eigenvalues is expected to be more involved due to 


the increased dimension of the system. 


5.2 FUTURE WORK 

Although the eigenvalue assignment techniques are con- 
venient from a computational viewpoint, there are serious dis- 
advantages for process control applications. This is mainly due to 
the fact that no systematic way of specifying the closed-loop eigen- 
vectors is available at the present time. The manipulation of 
available design freedoms in a trial and error manner is not always 
successful in a reasonable period of time and becomes more tedious 
as the number of states increase. 

In this respect, future work should be directed to the 
establishment of desirable closed-loop eigenvector configurations 
in state space. Since the closed-loop eigenvectors depend on the 
closed-loop eigenvalues and the design freedom is mainly provided 
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closed-loop eigenvectors, closed-loop eigenvalues and the set of 
vectors {g,} should be determined for the successful application 
of eigenvalue assignment techniques to the process control problem. 

Another importnat area where the future work is needed, 
is the minimization of gain elements in view of the possible limita- 
tions on the magnitude of available controls (i.e., control constraints). 
For multi-input systems, the gain elements are related to the selection 
of an arbitrary vector set, {g.}, and the relationship between the 
overall controller matrix and this set of vectors, {g,} would be 
a valuable asset for the control system design. 

These two unresolved problems in state feedback control 
are also important because the incomplete state or output feedback 
system design problem would profit from a fuller understanding of 
the state feedback problem. 

However, the simulation results should be verified 
experimentally before proceeding with a further investigation of 
eigenvalue assignment techniques. Unfortunately, the experimental 
evaluation of the simulated results for the evaporator system had 
to be postponed due to the equipment difficulties. Hopefully, the 


experimental studies will be possible in the near future. 
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NOMENCLATURE 


Alphabetic 


system matrix of a continuous-time model 


b control coefficient vector 

B control coefficient matrix 

BL first effect bottoms flow 

B, second effect bottoms flow 

C closed-loop system matrix 

Cc, first effect concentration 

C, product concentration 

Cy feed concentration 

d, disturbance element 

d disturbance vector 

D disturbance coefficient matrix or a desired closed-loop 
matrix defined in (3-30) 

EC ) ear polynomial of a square matrix 

f£ column vector defined in (2-3) 

g. element of vector g 

g column vector defined in (2-3) 

G square matrix defined in (3-33) 

F feed flow rate 

hy first effect enthalpy 

hy feed enthalpy 

H mode controllability matrix defined in (2-10) 

Hy mode controllability matrix defined in (3-17) 

H mode controllability matrix defined in (3-18) 
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i increment counter 

f identity matrix 

J Jordan canonical form of a square matrix 

k time counter 

K feedback control matrix or initial feedback control 


matrix defined in (3-11) 


proportional feedback control matrix 


=FB 

KF integral feedback control matrix 

m number of controls 

n number of states 

0 null matrix 

p number of eigenvalues to be changed 
q number of disturbances 

qg eigenvector updating coefficient 

Yr number of integrated states 

S steam flow rate 

iT discretization interval 

qT, first effect temperature 

Ty feed temperature 

I. coefficient matrix of integrated state vector 
u, control element 

u control vector 

u, control vector defined in (3-10) 

U, control vector defined in (3-10) 

Vv left eigenvector 


V left eigenvector matrix 


eo 
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right eigenvector 

right eigenvector matrix 
first effect holdup 
second effect holdup 
State variable 

State vector 

modal state vector 


integrated state vector 


element of vector OL 

column vector defined in (2-17) 
element of vector 6 

column vector defined in (2-17) 
discrete control coefficient matrix 
discrete disturbance coefficient matrix 
discrete state coefficient matrix 
open-loop eigenvalue 

set of eigenvalues 

eigenvalue matrix 

product 

desired eigenvalue 


sum 


time 
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Superscript 
. time derivative 
perturbation variable 


augmented vector or matrix 


7 recursive step counter 

— vector or matrix tanspose 

-] matrix inversion 

Subscript 

= vector 

7 matrix 

i element counter or run counter 
j element counter 

P size of matrix partition 
Symbol 

< > dot product 

{3 set 

Abbreviation 

det determinant 

diag diagonal 

P proportional control 

1a proportional plus integral control 


sign signum 


= 


he Se ea 


a 6 ay Py 


Ss 


+20% CF 


-20% CF 


+20F 


+207 TF 


steady state 

20% step increase in feed concentration 
20% step decrease in feed concentration 
20% step increase in feed flow rate 


20% step increase in feed temperature 
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APPENDIX 


DERIVATION OF MATRICES $(T), A(T) and 6(T) 


FROM THE CONTINUOUS AUGMENTED SYSTEM 
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If one directly discretizes the continouous augmented 
system given by Equation (3-2) following standard procedures in the 
literature [25], the coefficient matrices, $(T) , A(T) and 6(T) 


are expressed by: 


AT 
oGr) \=)¢e (A-1) 
2 ee 
A(T) = (f e dr)B (A-2) 
O 
x ek At a 
6(T) = (f e dt)D (A-3) 
oO 


I > 


where matrices, A, B and are defined by Equations (3-4), (3-5) 
and (3-6) respectively. If the matrix A in Equation (3-4) is non- 
singular, the right and left elgenyector matrices, W and V of 


matrix A can be related to the corresponding matrices’ W and V 


of matrix A by [32]; 


W = sag oe fecites editen (A-4) 
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F ae 
V = = en ma i -_- — = (A-5) 
ye aut paul 
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where A is a diagonal matrix whose diagonal elements are the eigen- 
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be expressed by the following relation, 


> 
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hi 
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em a Wey: (A-6) 
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- Lei ie’ 
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ee 


Substituting Equations (A-4), (A-5) and (A-7) into Equation (A-5) 


and rearranging the resulting expression gives, 


ES OW Hosta hen ees a ee meres (A-8) 


Hence, the integral expression in Equation (A-2) can be readily 


evaluated and given by, 


AT | 
-~_ = q 
a a A ae -1) BAC 
dt=|------f,---- - - '- - - (A-9) 
° eed 
ZA {A (e ~i)-Ti } L} II. 
From Equations (A-1) and (A-8), matrix $(T) is given by 
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After substituting Equation (A-9) into Equations (A-2) and (A-3) and 
using Equations (3-5) and (3-6), matrices ACT) and 8(T) are 


expressed by: 


_, 4 
4 A (eS 1 )8 
CON Gai teata eat Aieg moll (A-11) 
TA {A (e -I_)B-TB} 
AT 
: ae (en SLD 
OD aver one ee troll (A-12) 
an ae 
TA {A “(Ce -ZL)D-TD} 


From the equations given below Equation (2-33), it can be readily 


seen that: 


AT 

o(T) = e , (A-13) 
a7 ee 

ING ee Neb (A-14) 
=; 81 

Q(T) = A (e -L)D . (A-15) 


Substitution of Equations (A-13) - (A-15) into Equations (A-10) - 


(A-12) respectively gives the expressions given by Equations (3-37) - 


(3=39) 
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